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ABSTRACT 

Context. The so-called excess of cosmic ray (CR) positrons observed by the PAMELA satellite up to 100 GeV has 
led to many interpretation attempts, from standard astrophysics to a possible exotic contribution from dark matter 
annihilation or decay. The Fermi data subsequently obtained about CR electrons and positrons in the range 0.02-1 
TeV, and HESS data above 1 TeV have provided additional information about the leptonic content of local Galactic 
CRs. 

Aims. We analyse predictions of the CR lepton fluxes at the Earth of both secondary and primary origins, evaluate 
the theoretical uncertainties, and determine their level of consistency with respect to the available data. 
Methods. For propagation, we use a relativistic treatment of the energy losses for which we provide useful parameteri- 
zations. We compute the secondary components by improving on the method that we derived earlier for positrons. For 
primaries, we estimate the contributions from astrophysical sources (supernova remnants and pulsars) by considering 
all known local objects within 2 kpc and a smooth distribution beyond. 

Results. We flnd that the electron flux in the energy range 5-30 GeV is well reproduced by a smooth distant distribution 
of sources with index 7 ~ 2.3 — 2.4, while local sources dominate the flux at higher energy. For positrons, local pulsars 
have an important effect above 5-10 GeV. Uncertainties affecting the source modeling and propagation are degenerate 
and each translates into about one order of magnitude error in terms of local flux. The spectral shape at high energy 
is weakly correlated with the spectral indices of local sources, but more strongly with the hierarchy in their distance, 
age and power. Despite the large theoretical errors that we describe, our global and self-consistent analysis can explain 
all available data without over-tuning the parameters, and therefore without the need to consider any exotic physics. 
Conclusions. Though a standard paradigm of Galactic CRs is well established, our results show that we can hardly 
talk about any standard model of CR leptons, because of the very large theoretical uncertainties. Our analysis provides 
details about the impact of these uncertainties, thereby sketching a roadmap for future improvements. 
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1. Introduction 

I Cosmic ray (CR) electrons and positronfQ constitute ^ 1% 
. of the CR budget at the Earth in the GeV- TeV energy 
' range, and provide interesting means of probing the acceler- 
, ation processes in CR sources, propagation phenomenology, 
1 and the inter stellar environment itself, complementary to 
protons (e.g.. lBlandford fc Eichleijll987[ ). At energies > 100 
GeV, their observed properties are mostly set by the very 
local environment. Their typical propagation scale is indeed 
limited to the kpc scale because of the very efficient elec- 
tromagnetic energy losses caused by Compton scattering 
with the interstellar radiation fields (ISRF), the cos mic mi- 
crowa ve background (CMBV and the magnetic field (jJoned 
119651: iBlumenthal fc Gouldi Il970h . High energy CR elec- 
trons are produced directly by well-known astrophysical 
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CR accelerators such as supernova remnants (SNRs) and 
pulsars, in which case they are referred to as primary elec- 
trons. They can also be created by secondary processes, 
mostly nuclear interactions of cosmic protons and light nu- 
clei with the interstellar medium (ISM) gas concentrated 
in the Galactic disk (spallation) , in which case they are re- 
ferred to as secondary electrons. Because they have been 
assumed to hardly be produced in astrophysical sources, 
positrons have been proposed as potential tracers of new 
phy sics, in particular the a nnihilation or decay of dark mat- 
ter (|Silk fc Sredmckilll984[ ). Although the main theoretical 
ideas regarding the origin and propagation of cosmic elec- 
trons were forma lized a lon g time ago in the seminal mono- 
graph of .Ginzb m'g fc Svrovatsk ii (1964), their ongoing mea- 
surements are still far from being completely understood. 

T he observation by the PAM ELA satel- 
lite (jPAMELA Collaboration et all I2009D of a rising 
positron fraction up to ^100 GeV has triggered a consid- 
erable amount of interpretation attempts. Estimates of the 
cosmic electron and positron fluxes were first calculated 

^ Hereafter, electrons will denote both electrons and positrons, 
unless specified. 
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in detail in iMoskalenko fc Strond (jl998|) . where only 
secondaries were considered for positrons, that fail to 
match the PAMELA data. We derived novel predictions 
of the secondary positron flux at the Earth, with a 
particular focus on sizing the theoretical errors caused by 
uncertainties in spallation cross-sections, in the modeling 
of the progenitor interstellar (IS) CR flux, in the char- 
acterization of the energy losses, and in the propagation 
parameters (jPelahave et all l2009f ). Although the overall 
theoretical uncertainty is about one order of magnitude, we 
have still shown that a rising positron fraction was not ex- 
pected unless a very soft electron spectrum was considered. 
Even in that case, however, we have also illustrated how 
difficult it was to accommodate a good fit to the PAMELA 
data in spectral shape as well as in amplitude. This soft 
electron spectrum is at the lowest statistical edge of the 
current electron cosmic ray data below 30 GeV. Likewise, 
it is not supported by the unprecedented measurements 
performed with the Fermi satellite be tween 20 GeV and 1 
TeV of CR electrons plus positrons ()Abdo et al.ll2009bD . 
which sets the true denominator of the positron fraction. 
At this stage, separate data of positrons and electrons 
would be of particular interest and would provide stronger 
grounds to any interpretation attempt, but, unfortunately, 
are not yet available. Nevertheless, from both predictions 
of the secondary positron flux and the current data, it 
appears unlikely that this increase observed in the positron 
fraction is purely of secondary origin. Therefore, this 
positron excess points towards the existence of primary 
sources of positrons in the neighborhood. Note finally 
that the cut-off in the electron flux observed by HESS 
around 3 T eV provides interestin g and complementary 
information (jAharonian et al.]|2008l ). 

It has long been demonstrated that astrophysical 
sources may supply this ex tra-yield of cosm ic positrons. 
For instance, as d iscuss e d by Boulares (I1989D fsee also e.g. 



Aharonian et al.l Il995l : IChi et a". 

2001[) , pulsars could provide sizable contributions to 
positron flux from pair conversions of 7-rays in 



19961 : IZhang fc Chend 
the 
the 



conversions 

strong magnetic fields that they host. 
cently revisited by s everal authors ( e.g. iHooper et al.ll200' 



This has been re- 



lYiiksel et all 120091: iProfumol l2008l : iMalvshev et al.l 1200' 
who have drawn similar conclusions. Another class of con- 
tributions invokes spallation processes with the ISM gas 
during the acceleration stage of cos mic rays inside SNRs 
that had n o t been considered before (IBerezhko et al.ll2003 ; 



Blasill2009l: felasi fc Serpicol 120091: IMertsch fc Sar^ 



2006; 



Ahlers et al. I I2009D . This hypothesis leads to the addi 



tional production of antiprotons or heavier secondary nu- 
clei, providing interesting counterparts that should be ob- 
served in the near future. Finally, using a more refined spa- 
tial distribution of sources and interstellar gas might also 
lead t o a rising positron fraction in the PAMELA energy 
range ijShaviv et al.ll2009t) . 

The large amount of standard, but still different, as- 
trophysical interpretations of the observed positron frac- 
tion is noteworthy and points chiefly towards significant 
lacks in our understanding of the cosmic electron pro- 
duction in sources and their subsequent propagation in 
the Galaxy. This also points towards a standard model of 
Galactic cosmic rays still being far from complete, in spite 
of the progresses achieved so far in the description of cos- 
mic ray sources, propagation and interaction with the ISM 
and ISRF. Because high energy electrons have a propa- 



gation horizon much smaller than light cosmic ray nuclei, 
they offer an interesting means to improve the overall phe- 
nomcnological modeling, the Galactic environment being 
indeed much more tightly constrained locally. We never- 
theless emphasize the robustness of the standard paradigm 
of cosmic ray physics as formali zed in the seminal book 
of iGinzburg &: Svrovatskiil (jl964 ): distinguishing standard 
model from standard paradigm appears to us important to 
avoid considering departures from peculiar observational 
data as deep failures of generic astrophysical explanations. 

Our purpose is to develop novel calculations of the elec- 
tron and positron fluxes to assess the relative roles of the 
different primaries and secondaries in the positron fraction. 
This is somehow a continuation of the study that we per- 
formed on secondary positrons (jPelahave et al.l [20091 ). We 
treat all of these components in a self-consistent framework 
that includes e.g. improved propagation modeling (with 
a full relativistic treatment of the energy losses) as well 
as constrained properties of local sources, including both 
SNRs and pulsars. In addition to improving and clarify- 
ing the interpretation of the PAMELA data from standard 
astrophysical processes, this study also helps us to verify 
whether the cosmic positron spectrum can provide inter- 
esting perspectives in the search for new physics. A partic- 
ularly important issue is whether positrons injected from 
dark matter annihilation could be differentiated from all 
other astrophysical contributions. Dark matter could in- 
decd in some cases manifests itself in this channel (e.g. 
iBahz fc Edsial 119981: iHooper fc KriQ 120041 iLavalle et al. 



2007':'Asano et al.ll200 7^'Bcrgst rom et al .'200g:'CireUi ct al 



2008; Delahave et al.l 12008; Fieri et al.l 12009: Catena ct al. 
l2009l ). and the discovery of an exotic contribution to the 
positron budget would be a spectacular result. Any such 
result would, however, have to rely on solid grounds, in 
particular a good understanding of the astrophysical con- 
tributions. We show that the theoretical uncertainties are 
very large, and discuss in detail the relative impact of each 
ingredient. This variance in the predictions is quite bad 
news for exotic searches, because it indicates that the back- 
ground is poorly constrained. Despite these uncertainties, 
we show that our calculations, involving pure astrophysi- 
cal processes in a self-consistent framework, can very well 
explain the whole set of available data on CR leptons, with- 
out over-tuning the parameters, and therefore without any 
need of exotic physics. 

The outline of the paper is the following. In Sect. |5J we 
describe in detail our propagation model, with a particular 
focus on the relativistic treatment of the energy losses. In 
Sect. 121 we revisit the predictions of the local secondary 
electron and positron fluxes and discuss the effects of our 
improved propagation model co mpared to the results we 
derived in iDelahave et al.l (|2009t) . In Sect. (H we compute 
the primary electron component by considering a smooth 
distribution of SNRs beyond ~ 2 kpc from the Earth, and 
by determining the contribution of each known SNR within 
this distance; we also discuss in detail the impact of the 
source modeling. In Sect. [SJ we briefly revisit the primary 
positrons that pulsars could generate by using the same 
approach as for electrons. We flnally compare our results 
with all available data on CR leptons in Sect. |S1 before 
concluding in Sect. [71 
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2. Propagation of electrons and positrons 

2.1. General aspects 

CR propagation in the Galaxy involves quite complex 
processes. The spatial diffusion is caused by convection 
upwards and downwards from the Galactic disk and by 
the erratic bouncing of CRs off moving magnetic inho- 
mogeneities, which also induces a diffusion in momentum 
space, more precisely diffusive reacceleration. Energy losses 
along the CR journey have additional effects on the diffu- 
sion in momentum space. Nuclei can also experience nu- 
clear interactions (spallation) ; this is of course irrelevant to 
electron propagation, but spallation will still be considered 
as the source of secondaries. The propagation zone spreads 
beyond the disk, and is very often modeled as a cylindrical 
slab of radius R ~ -Rdisk — 20 kpc, and a vertical half-height 
of L ~ 1 — 15 kpc. Astrophysical sources of CRs and the 
ISM gas are mostly located within the disk, which has a 
vertical extent of Zdisk — 0.1 kpc. More d etails on propaga- 
tion p henomenology can b e found in e. q. Bcrcz inskii et ahl 
(ll990D .rLongair (199j) and lStrong ei~al. (2007). 

Throughout this paper, we discuss high energy electrons 
particularly of energies above ~ 10 GeV, for which the 
effects of solar modulation are much weaker. We demon- 
strated in iDelahave et al.l ( 2009f ) that convection and reac- 
celeration could be neglected above a few GcV, so that 
the propagation of electrons can be expressed in terms of 
the usual current conservation equation VM = Q{E,x,t), 
where the transport operator V can be expanded as 



rlK 

dtN-V- {K(E)VN} + dE { ^Af 



Q{E,x,t). (1) 



The electron number density per unit of energy is denoted 
Af = Af{E,x,t) = dn/dE, K{E) is the energy-dependent 
diffusion coefficient assumed isotropic and homogeneous, 
dE/dt is the energy- loss term and Q is the source term. As 
mentioned above, we neglected convection and reaccelera- 
tion. 

The above equation can be solve d numerically, e.g. by 
mean s of the public code GALPROP (IStrong fc Moskalenkol 
|1998() . which treats CR nuclei and electrons in the same 
framework. However, most of the studies using this code 
do not usually correlate the f eatures of protons at sources 
with those of elec trons {e.g. iMoskalenko fc Strong! Il998t 
IStrong et al.ll2000l ). which alleviates the relevance of treat- 
ing nuclei and electrons in the same global numerical frame- 
work. In that case, one can always tune the source modeling 
differently for each of these species to accommodate the ob- 
servational constraints. 

We adopted instead a semi-analytical propagation mod- 
eling, which is designed to survey a wider parameter 
space and both clarifies and simplifies the discussion 
on theoretical uncertainties. Analytical steady-state solu- 
tions to Eq. (TTI). in terms of Green f unctions, can be 
found in e.g. 'Bulanov & Dogcl' (ll974^. iBerezinskii et al.l 
diggol). iBa hz fc Edsio (1998). iLavalle et al.l (12001 7 
or Delahave et al.l (j2008f ). in the non-relativistic Thomson 
approximation of the inverse Compton energy losses. We 
improve this model by including a full relativistic calcu- 
lation of the energy losses (see Sect. 12. 4p and the time- 
dependent solution to Eq. ([ij, which has to be used when 
dealing with local sources (see Sect. 12. 2|) . The propagation 
parameters are constrained as usual, by means of the ra- 



tio of secondary to primary stable nuclei, except for the 
energy-loss parameters, which are constrained from the de- 
scription of the local ISRF and magnetic field. This latter 
point is discussed in detail in Sect. 12.51 

For the sake of completeness, we briefly recall the Green 
functions that are steady-state solutions to Eq. ([T]), disre- 
garding the energy-loss features for the moment. Assuming 
that spatial diffusion and energy losses are isotropic and ho- 
mogeneous, it is an academic exercise to derive the steady- 
state Green function in an infinite 3D space, which obeys 
Vtg^S^{x-Xs)S{E~Es), I.e., 



g{x,E^Xs,E,) = 



1 



b{E) (7rA2) 



r-expi^^^i , (2) 



where the subscript s flags quantities at source {Eg > E), 
and we define the energy-loss rate and the diffusion scale 
to be 



, dE 



,,KiE') 



4 / dE',, 



(3) 



The propagation scale A characterizes the CR electron hori- 
zon and depends on energy in terms of the ratio of the dif- 
fusion coefficient to the energy-loss rate. If these are both 
described by power laws, e.g., K{E) cx E^ and b{E) oc i?", 
then g cx i^f "iC'^+i); this is of importance when discussing 
the primary and secondary contributions later on. For def- 
initeness, we define 



K{E) 
KE) 



/3A'o 



iGvy 

Eo 



Koe' 



E 



e" with e = 

Ti Ea = l GeV 



(4) 



where Kq and bg are the normalizations of the diffusion co- 
efficient and the energy-loss rate, respectively, that carry 
the appropriate dimensions, and r; is the characteristic 
energy-loss time. 

Because of the finite spatial extent of the diffusion slab, 
boundary conditions must be taken into account when the 
propagation scale is on the order of the vertical or ra- 
dial boundaries. At the Earth location, which we fix to 
be {xQ,yQ, zq) = (8,0,0) kpc throughout the paper, the 
radial boundary is irrelevant while (i? — tq) > L, which 
is almost always the case for reaso nable values of L and 
R, constrained by observ ations (e.g. IStrong fc Moskalenkol 
119981: iMaurin et al.ll200ll ). Therefore, we briefly review the 
solutions accounting for the vertical boundary condition 
only. In that case, one can split the general Green function 
into two terms, one radial and the other vertical, such as 
G = {Gr X Gz)/b{E). The radial term is merely the infinite 
2D solution 



Grir,E<~r,,Es) = ^cxpj- ^'' J'^' 



(5) 



where r is the projection of the electron position in the 
z = plane, and the subscript s refers to the source. The 
vertical solution can be determined by different methods. 
On small propagation scales, more precisely for A < L , 
one can use the image method {e.g. ICowsik fc Lei 119791: 
iBaltz fc Eds"ial[T998h 



Gz{z,E^ Zs,Es) 



+ 00 

E 

n— — oo 



(-1)" 
^/^A 



GXp 



A2 



,(6) 
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where Zs^n = 2nL + ( — l)"zs. On larger propaga- 
tion scales, i.e. A > L, the basis defined by the 
Hclmholt z eigen-functions a llows a better numerical con- 
vergence (iLavalle et al.l[2007l ). In that case, we have instead 



n=l (. 



(7) 



where the pair and odd cigen-modcs and eigen-functions 
read, respectively 

kn — (n — l/2)n/L ; fc^ = mr / L ; 
0„(z) = sin(fc„(i - |z|)) ; (j)'^{z) = sin {k'^{L - zj) . (8) 

The radial boundary condition becomes relevant when {R— 
tq) ^ L ^ X. We accounted for it using the image method 
for the radial component, or, since the smooth source term 
exhibits a cylindrical symmetry, b y expanding the solution 
in terms of Bessel series (see e.g. Bulanov fc Dogel 119741 : 
iBerezinskii et alj 119901: iDelahave et alj I2008D . The radial 
boundary condition is, however, mostly irrelevant in the 
following, since we will mainly consider electron energies 
> 10 GcV, for which the propagation scale is no more than 
a few kpc. 

2.2. Time-dependent solution 

The steady-state solutions derived above are safe approxi- 
mations for a continuous injection of CRs in the ISM, as in 
the case of secondaries. In opposition, primary CRs arc re- 
leased after violent and localized events such as supernova 
explosions, the remnants and sometimes pulsars of which 
are assumed to be the most common Galactic CR acceler- 
ators. Since the supernova explosion rate is most likely 
a few per century, the CR injection rate could exhibit sig- 
nificant local variations over the CR lifetime (confinement 
time, or energy-loss time, depending on the species) pro- 
vided this latter is much longer than the individual source 
lifetime. Since electrons lose energy very efficiently, there is 
a spatial scale (an energy scale, equivalently), below (above) 
which these local variations will have a significant effect on 
the local electron density. To roughly estimate this scale, 
one can compare the energy loss rate b{E) with the local 
injection rate. Assuming that source events are all identical 
and homogeneously distributed in an infinitely thin disk of 
radius i? = 20 kpc, local fiuctuations are expected to be 
smoothed when integrated over an electron horizon A such 
that r^(A/i?)2 > b{E)/E. Using Kq ~ O.OlkpcVMyr, 
b{E) w (GeV/Myr/315)e2 and w 1/100 yr, we find 
that E 80 GeV, which means that local fluctuations 
of the flux are probably important above a few tens of 
Ge V. A s i milar reasoning was presented a few decades ago 
by IShenI (Il970h . We recall that a significant number of 
SNRs and pulsars are actually observed within a few kpc 
of the Earth. Therefore, current multi-wavelength measure- 
ments may help us to feature them as electron sources, and 
thereby predict the local electron density. 

To estimate the contribution of local transient sources, 
we need to solve the full time-dependent transport equa- 
tion given in Eq. ([T]), and we demonstrate that the method 
used for the steady-state case can also be used, though 



partly, for the transient case. The time-dependent Green 
function, Qt, is defined in terms of the transport oper- 
ator, asking that VQt = 5^{x — Xs)S{E — Es)S{t — ts). 
To so lve this equation, we generally work in Fourier space 
(e.g. iGinzburg fc Svro vatskiill l964: IBere zinskii et al.lll990l: 
Atoyan et al.l 119951: iKobavashi et al.l 12004 : Baltz fc Wail 



20041 ) ■ using 



Qt{t,E,x) 



d^kduj 

X exp {i(fe • a; -f tot) } 0(a;, E, fc) 



(9) 



In Fourier space, we derive the ordinary differential equa- 
tion for E, for each pair {uj, fc) 



{iu + k^K{E)} (j)-dE (b{E)4>^ = S(E - E,) 



(10) 



1 



(27r)2 

which is solved by the function 

0(a;,i?,fc)=^^exp|-ifc^A2 
X exp{-'iijj{ts + Ar)} . 



exp {—i{k ■ Xs + wis)} 



%k ' X t 



(11) 



This solution is only valid for E < Eg because it describes 
processes ruled by energy losses. It contains the propaga- 
tion scale A previously defined in Eq. ([3]) and the loss time 
defined as 



At{E,E,) 



dE' 
E b{E') 



(12) 



This loss time corresponds to the average time during which 
the energy of a particle decreases from Eg to E because 
of losses. The inverse Fourier transformation is straightfor- 
ward from Eq. and we eventually obtain 



gt(t,E,x ^ ts,Es,Xg) 



5{At - At) 



HE) 



(7rA2)3/2 



(13) 



where At = t — tg. We recognize the product of the steady- 
state solution and a delta function mixing real time and loss 
time. As in the steady-state case, we can account for the 
vertical boundary condition by expanding this 3D solution 
by means of the image method or on the basis of Hclmholtz 
eigen-functions. The final result can therefore be expressed 
in terms of the full steady-state solution 

gtit,E,x ^ tg,Eg,Xg) ^ SiAt- AT)g{E,x ^ Eg,Xg) . 

(14) 

An alternative interpretation of the time dependence comes 
up when the temporal delta function is converted into an 
energy delta function, which has been shown to be appro- 
priate for bursting sources for which At is fixed. In this 
case, the Green function is instead given by 



gt{t,E,x^ts,Es,Xg) = S{Es^E*)b{E*) 

xg{E,x^ Es,Xs 

where the energy E* satisfies 

At{E,E*) = At. 



(15) 
(16) 
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Thus, E* corresponds to the injection energy needed to ob- 
serve a particle with energy E after a time At = t — ts- 
Although there is no analytical solution to this equation 
in the full rclativistic treatment of the energy losses (see 
Sect. I2.4p . we can still derive it in the Thomson approxi- 
mation 



E-" 



with EZ^^ 



Th. 



E 



approx. 1-E/EZl 



(17) 



[boAtr'^^^Eo 



where we used the energy- loss term from Eq. (|4]). We see 
that when the energy-loss timcscalc r; ^ At, we have E* ss 
E. Wc also see that a maximal energy is set by the ratio 
Ti/At: in the Thomson approximation, a particle injected 
with energy > E^^^ will have already lost all its energy 
by At. We emphasize that i?niax > -E'max ™ the general 
relativistic case (see Sect. 12. 4|) . 

Note that an additional consequence of this energy E* 
is that the propagation scale A is no longer set by energy 
losses but instead by the injection time At. In the simpli- 
fied case of a constant diffusion coefficient K, we would 
indeed have found that = AK At. Of course, the energy 
dependence of the diffusion coefficient slightly modifies this 
relation, but this remark will further help us to make a 
rough prediction about the observed spectrum for a burst- 
ing source (see Sect. 12. 3p . 

Finally, we underline that solutions to the time- 
dependent transport equation do not always provide causal- 
ity, which is important to avoid incorrect predictions when 
varying the source age and distance. To ensure causality at 
zeroth order and for the sake of definitivcness, we use 

gt{t,E,x^ts,Es,Xs) - e{cAt-\\x-Xs\\)6{Es- E*) 

xb{E-')g{E,x^ E„Xs) , (18) 

as our time-dependent propagator. A more accurate 
causal solution would need more specific methods inferred 
from e.g. detailed studies of relativistic heat conduction. 

2.3. Approximated links between propagation models and 
observed spectra 

To anticipate the discussion about the observed versus pre- 
dicted spectra primary and secondary electrons, it is useful 
to show how observed indices can be formally linked to the 
propagation ingredients. We now establish approximate re- 
lations between the observed spectral index 7, the source 
spectral index 7, and the propagation parameters. In the 
most general case, the interstellar (IS) flux at the Earth, 
i.e. without accoimting for solar modulation, is expressed 



47r 



(19) 



dtsdEsd^Xsg{E,XQ ^ ts,Es,Xs) Q{ts,Es,Xs) 



We first discuss the steady-state case, before continuing to 
the case of time-dependent sources. 

The energy dependence arising in the electron propaga- 
tor (see Sect. 12. ip comes from spatial diffusion and energy 
losses. At high energy, one can assume that the propaga- 
tion scale is short enough to allow us to neglect the vertical 



boundary condition, so that one can use the 3D propaga- 
tor to predict the electron flux on Earth, given a source 
Q(E,x). Since we consider a short propagation scale, and 
since sources are located in the Galactic disk, we can as- 
sume a source term that is homogeneously distributed in 
the disk. This is a very good approximation for secondaries 
(see Sect. [3]), and fair enough for primaries (see Sect. 14. 2|) . 
Likewise, we consider that the source spectrum is a mere 
power law of index 7, so that the source term can be writ- 
ten as Q{E, x) = 2hQo S{z) €~'^, where h is the half- height 
of the disk and e, which is deflned in Eq. is the dimen- 
sionless energy parameter. Given this source term, the flux 
on Earth is given by 



ME) 



(20) 



where, o = \/a — S — l/{-f — 1) = 0(1), and we have used 
the 3D propagator defined in Eq. ([5]), the energy depen- 
dence of which is fully determined from Eqs. ^ and (j?]). 
Accordingly, the spectral index 7 after propagation reads 



1^1+ -{a + 5-1) 



(21) 



As discussed in Sect. 12.41 the energy-loss rate is dominated 
by inverse Compton and synchrotron processes. In the non- 
relativistic Thomson approximation, we have a — 2, lead- 
ing to7 = 7-|-^(l-|-5). From this basic calculation, it is 
easy to derive rough values for 7 and 6 consistent with any 
spectral index 7 measured on Earth. For instance, 7 « 3 
translates into a source index 7 in the range [2.1,2.35] for 
S e [0.3,0.8]. Although very useful to first order, this crude 
spectral analysis is only valid for a smooth and flat dis- 
tribution of source(s), and signiflcantly differs when local 
discrete effects are taken into consideration. Implementing 
full relativistic losses induces a = agg ^ 2, which implies a 
harder 7. This will be delved into in more detail in Sect. 21 

Finally, we extract the observed spectral index 7^ for 
a single event-like source, which differs slightly from the 
above calculation. The source term can be expressed as 

= Q*,oS{\xs \ — d) S{ts — t*) — we discuss this injec- 
tion spectrum in more detail in Sect. 14.31 Assuming further 
that the source is located within the propagation horizon 
d ^ X and that a burst occurs at a time much earlier than 
the energy-loss timescale t^, -^ti, we readily flnd that 



, Pcb{E^) Q.,oe*^ 



47r h{E) (7rA2)3/2 47r (4 tt t^)3/2 



where 

1. 



3. 



(22) 



(23) 



Here, we have considered that the propagation scale is no 
longer fixed by energy losses, since ^ r;, but instead by 
(see the discussion at the end of Sect. 12. 2p . In this case, 
since i?* ~ E, the spectral index is not directly affected by 
the energy losses. 

2.4. Full relativistic energy losses 

In the GeV-TeV energy range, electrons lose their en- 
ergy by electromagnetic interactions with the ISRF (inverse 
Compton scattering) and the magnetic field (synchrotron 
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emission), while Bremsstrahlung, ionization and Coulomb 
interactions with the ISM are negligible. Most studies have 
used the Thomson approximation to account for inverse 
Compton losses, which is valid for an electron Lorentz fac- 
tor 7e S, TTieC^ / E^\i, whcrc - E ph is the photon energ y {e.g. 
iMoskalenko fc Strongj Il998t fPelahave etlll 120091 ). This 
translates into a maximal electron energy of ~ 1.11 x 10® 
GeV for interactions with CMB (E'ph ~ 2.35 x 10"'^ cV), 
and of 7.58 x lO-^/S.ee x 10^ GcV for IR / starlight radi- 
ation, respectively (with -EphjR/* — 3.45 x 10~'^/0.3 cV). 
From those numbers, it is clear that the Thomson approxi- 
mation is no longer valid for energies at Earth above a few 
tens of GeV, for which a full relativistic description of the 
term dE/dt in Eq. ([T]) is consequently necessary. Few other 
studies have imp lemente d this relativistic treatment (e.g. 
iKobavashi erall[2004 : Schhckeiser fc Ruddc I 20091. 

The calculation of inverse Compton scattering of elec- 
trons with photons in the relativistic regim e was de- 
rived in the astrophysical context by iJoned (|l965f ). It 
was subsequentl y extens i vely r evisited and complemented 
by [glumenthal fc Gouldl ()1970l ) . In the following, we rely 
on the latter reference to derive our relativistic version of 
the inverse Compton energy losses, to which we refer the 
reader for more details. 

We consider relativistic electrons propagating in an 
isotropic and homogeneous gas of photons, which, more- 
over, exhibits a black-body energy distribution. The rele- 
vance of these assumptions will be discussed in Sect. 12.51 
The electron energy-loss rate can be expressed in terms of 
the energies e and ei of a photon before and after the col- 
lision, respectively, as 



dE 
~dt 



= / dt 



j dei(ei 



dN, 



coll 



dt de dei 



(24) 



The collision rate is given by 
3 (Tt c dn{e)/de 



dNr. 



dt de dei 



47e e 
2q{\nq-q 



2 (1 + Tq) 



(25) 



where dn{e)/de is the initial photon density in the energy 
range de, which, for black-body radiation has the form (in- 
cluding the two polarization states) 



dn 



47re^ 



— 2 X 
de {2Tihcf 

and 



£1 



ei 



r = 



47ee 



(26) 



(27) 



r(l - ei) 7emc" 

From kinematics, the range for ei is readily found to be 
[e, |Y^py], which translates into [-^, 1] for q. It is convenient 
to rewrite the energy loss rate in terms of an integral over 



dE_ 

dt 



de I dq 



r2(7emc2)2 ( q 1 
(l + rg)2 \ (1 + Fg) " 4^ 



dt de dei 



(28) 



where the integral over g is found to be analytical, so that 
one can easily check the full numerical calculation. 



We define a dimensionless parameter that characterizes 
the relevant regime to be used for the energy loss rate 



7e (fefcTb) 



(29) 



where Tq is the mean temperature of the radiation field. 

The non-relativistic Thomson limit is recovered for in- 
verse Compton processes within a black-body radiation 
field, using F <C 1 or equivalently a ^ 1 



dE 
'dt 



-CTTCf/radTe 



(30) 



where ?7rad = / deedn/de, whereas the Klein-Nishina 
regime applies for a ^ I 



dE {rrii^ckhTo)'^ 
'dt 



16 



In 



47e fcfcTo 



1.9805 



(31) 



In Fig.[l] we compare both regimes with the full calculation. 
From our numerical results, we derived a parameterization 
that is valid for any black-body radiation field, given by 



{Thomson for Cn-r 

^^exp{E..o^.(ln«)^} forC., (32) 
Klein — Nishina 



-'int 
for Cu-r 



where the conditions C read 



r 

Cint 

r 



a < 3.8 X lO^'' , 

3.8 X lO"'' < a < 1.8 X 10^ 

a > 1.8 X 10^ . 



(33) 



The fitting formula associated with the intermediate regime 
provided in Eq. p2p may be used with the parameters 



Q = |74.77, -0.1953, -9.97 X 10~^ 

4.352 X 10"^, 3.546 x lO"'^, -3.01 x 10"^ 



} ■ (34) 



An additional smooth interpolation between these three 
regimes might improve the calculation by avoiding tiny gaps 
at connections, which could arise e.g. from very small nu- 
merical differences in the unit conversions or constants used 
above. This parameterization is valid for any black-body 
distribution of photons. If one considers a black-body dis- 
tribution, the absolute energy density J7iad of which differs 
from that generally derived U^^^ , then one can simply renor- 
malize Eq. ((32|) by a factor Urad/U^ad ^'^ obtain the correct 
energy-loss rate. 

In the following, we use Eq. ([5^ to describe the energy 
loss rates associated with Compton processes. 



2.5. Review of the propagation parameters 

In addition to energy losses, five parameters regulate the 
diffusion properties of Galactic CRs: Kq and S defining 
the diffusion coefficient (see Eq. 0]), the half-thickness of 
the diffusion zone L, the convective wind velocity Vc, and 
the Alfven speed of magnetic field inhomo geneities Vg, re- 
sponsible for reacceleration. It was shown in lDelahave et all 
(120I39I ) that the last two effects can be neglected above a few 
G eV. These paramete rs were self-consistently constrained 
in lMaurin et al.l (j200lt ) with ratios of secondary to primary 
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a s y T / ITlg 



Fig. 1. Comparison between the different relevant regimes of 
the inverse Compton energy loss for any black-body radiation 
field. 



Model 


S 


Ko 


L 




Va 






[kpcVMyr] 


[kpc] 


[km/s] 


[km/s] 


min 


0.85 


0.0016 


1 


13.5 


22.4 


med 


0.70 


0.0112 


4 


12 


52.9 


max 


0.46 


0.0765 


15 


5 


117.6 



Table 1 . Beacon sets of diffusion parameters derived in 
iDonato et al.l (120041') compatib le with the B/C analysis per- 
formed in iMaurin et al.l (|200lD . The med setup will be our de- 
fault model. 



nucle i — mostly boron to carbon B/C (see also lPutze et al.l 
1201(1 for a more recent analysis). In the following, we use 
the available parameter space provided by these authors. 
Nevertheless, as useful beacons for brackctting the theo- 
retical uncertainties, we also use the mm, med, and max 
subsets of propagat ion parameters, which were derived in 
iDonato et al.l (j2004 ) and called so after the hierarchy found 
on the primary antiproton fluxes, for sources spread all over 
the diffusion zone (not only in the disk). These models arc 
described in Table [T] 

For the normalization of the diffusion coefficient Kq, 
we note that B/C measurements actually constrain Kq/L, 
not Kq alone (jMaurin et al.ll200ll) . Moreover, using radioac- 
tive species does not all ow yet to fully break this degener- 
acy (IDonato et al.l[200^ . This explains why the min (max) 
configuration, which has a small (large) L, is associated 
with a small (large) value of Kq. The spectral index of the 
diffusion coefficient S decreases from min to max, which is 
important to the spectral analysis of the electron flux. 

Tighter constraints are expected to be possib le with fu- 
ture PAMELA data, and hopefully with AMS2 (jBattistonl 
|2007() . The current uncertainty in those parameters 
leads to large theore tical errors in secondary positrons 
(jPelahave et al.l [20091 ) and therefore electrons, as reviewed 



in Sect. [31 The error in astrophysical primaries is assessed 
in Sect. 1131 

In contrast to stable nuclei in the GcV-TcV energy 
range, electrons are strongly affected by energy losses, 
which ha ve a significant effect on their transport. It was 
shown in iDelahave et al.l ()2009() that inverse Compton and 
synchrotron processes dominate in this energy domain. 
Therefore, it is crucial to constrain the ISRF — includ- 
ing the CMB, and radiation from dust and stars — and the 
magnetic field as accurately as possible, within the horizon 
of GeV-TeV electrons, i.e. 0(1 kpc). 

In Sect. 12.41 we developed a method to calculate the 
inverse Compton energy losses in a fully relativistic formal- 
ism, provided the target radiation fields can be described 
in terms of black-body distributions. This is obviously 
the case for the CM B, the tempe rature of which was re- 
cently re-estimated in lFixsenI (|2009D to be 2.7260±0.0013 K. 
However, it is well known that the ISRF is not simply 
Planckian radiation, since it consists of many different com- 
ponents — IR radiation from dust, optical and UV radia- 
tion from stars, diffuse X-ray emission, etc. — with different 
spatial distributions. Since by using CMB only we estimate 
the electron propagation scale to be < 2 kpc for electron en- 
ergies above 10 GeV, we disregard the spatial dependence 
of the ISRF, and only consider local averages. 

In the left panel of Fig. |21 we re port the ISRF dat a 
that we extracted from the analysis of iPorter et al.l i 



and averaged in cylinders of radius and half-height of 500 pc 
(model 2, M2) and 2 kpc (model 1, Ml) about the Earth, on 
top of which we show that a sum of black-body distributions 
can provide a reasonable fit. These two models, defined with 
a set of components characterized by their temperatures 
and energy densities, are summarized in Table O They can 
be used to estimate the theoretical error coming from uncer- 
tainties in the characterization of the ISRF. We may assume 
that these uncertainties reflect those affecting the data that 
we used, though error bars are not available. Note that the 
parameterizations appearing in Tabic ID arc not designed 
to reflect the true radiative physics operating in the ISM, 
which is beyond the scope of this paper. Nevertheless we 
observe an interplay between the IR and UV components, 
depending on the averaging volume: taking a smaller vol- 
ume results in a larger (smaller) IR (UV) contribution due 
to the efficient UV-absorption and IR-cmission properties 
of the dust, which is mostly concentrated in the disk. 

The synchrotron emission can also be expressed as an 
inverse Compton scattering on a black-body distribution of 
virtual photons from the magnetic field. In this case, the 
characteristic energy of the radiation field is given by the 
cyclotron frequency 



h Ur 



heB 
2 TT m 



1.16 X 10 



-14 



B 



1/iG 



cV , (35) 



where m is the electron mass, ks is Boltzmann's constant, 
and B is the value of the Galactic magnetic field. It is 
clear that the condition 7 Eb ^ rnc^ is fulfilled for the 
whole electron energy range considered in this paper, so 
that the Thomson approximation is fully valid. We esti- 
mate the local magnetic field relevant to th e synchrotron 
losses to i?sync ~ l^^G (scc e.g. lFerrierel200lD . for which the 
corresponding energy density derived from classical electro- 
dynamics is Ub = Bsync/(2/^o)- 

The synchrotron energy losses do not depend on the 
mean value of the magnetic field (B), but on the mean 
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Thomson reg. aprox. 

Model 1 
Model 2 
CMB 

CMB + Synchi 



Model 2 - 
CMB - 

CMB + Synchrotron - 



10 



10' 



10' 




Fig. 2. Left: Energy densit y distribution of the ISRF averaged in 2 boxes of different volumes, where all components appear (data 
taken from the analysis by iPorter et alj (|2008l ')'). Our models 1 and 2, using both black-bodies for all components, are reported 
against the data. Middle: Corresponding energy loss rate. Right: Toy electron fluxes associated with the previous energy loss 
configurations, assuming an injection spectrum oc E^^. 







To[K] 


U,,d [10-'' GeV.cm--^] 




CMB 


2.726 


Planckian (b-b) 


1—1 


IR 


33.07 


25.4 (4.5 10"'' X b-b) 


Stellar 


313.32 


5.47 (1.2 lO^^x b-b) 


UV 


3,249.3 
6,150.4 
23,209.0 


37 (7.0310-"x b-b) 
22.9 (3.39 10"" X b-b) 
11.89 (8.6710""x b-b) 


CM 


IR 


33.653 


32.12 (5.310"''x b-b) 


Stellar 


313.32 


6.2 (1.36 lO-^'x b-b) 


UV 


2,901.13 
5,570.1 
22,048.56 


33.76 (1.01 10"'^ X b-b) 
25.93 (5.7 10"" X b-b) 
10.26 (9.210-''^x b-b) 



Table 2. Parameters used to fit the local ISRF with black- 
bodies. They corres pond to fits perform ed on the data extracted 
from the analysis of lPorter et al.l (|2008l l and averaged over cylin- 



ders of radius and half-height of 2 kpc (model Ml) and 0.5 kpc 
(model M2) about the Earth. We choose Ml as our defauh ISRF 
model. In the last column, (ax b-b) indicates that the component 
can be obtained from a standard black-body spectrum renormal- 
ized by a factor of a. 



value of the squared field (S^). Although the mean value 
(B) = Br, namely the regular component of the magnetic 
field, the irregular component implies that (_B^) > {B) . 

iJafFe et all (|201Cl[ ) provide constraints on the different 
components of the magnetic field inside the disk. There arc 
actually three components: the regular component Br and 
two irregular components (one aligned with the regular one 
Ba, the other completely isotropic Bi). The relevant value 
we need for the synchrotron losses is 



Ba^l-ir 



'Br 



(36) 



where it is easy to show that the variance is given by cr^ 
{Bl) + {Bf). (Note that (S^) = {Brf). 



From the results obtained bv lJaffe et al.l (|201Cl[ ). we find 

(37) 



that 

(Br) = 1-3/iG, 

Bf) = 2.1 - 4.2 /iG [intcrarm - ridge] , 



B^„ 



0.0 



3.3 fiG [intcrarm - 
which translates into the range 
2.32 < Bsync < 6.13 AiG. 



ridge] 



(38) 



Nevertheless, we should not neglect the vertical dependence 
of the magnetic field, which is usually found to be expo- 
nential, with a typical scale of ~ 1 kpc. We emphasize that 
this scale is generally obtained from Faraday rotation mea- 
sures as well as from observations of the all-sky polarized 
synchrotron emission, which translate into magnetic field 
intensity only after deconvolution of the thermal and non- 
thermal electron density. This density is usually grossly 
modeled by assuming the sum of Boltzmann and single- 
index power-law spectra, and a z-cxponcntial spatial de- 
pendence, which is itself motivated fro m radio observations 
of the same synchrotron emission (see iJansson et all 120091 
for a recent analysis). Therefore, these estimates might be 
affected by potentially large systematic errors. 

We assume for simplicity that all components are con- 
stant in the disk (justified at the kpc scale around the 
Earth) and have the same vertical behavior, so that 



-^sync ^) 



Soexp<i 



If we average Sgync inside a spherical volume of radius 2 
kpc, which corresponds to the typical propagation scale for 

electrons, we find that ^J'{B^^^ = 0.05 (0.37) x Bq assum- 
ing that Zq = 0.1 (1.0) kpc. Given the range in Eq. ([55)1 . we 
obtain 



0.11 mG< (Ssync) 
0.87 mG< (Ssync) 



< 0.30 /xG 

< 2.29 /xG 



[zq 



0.1 kpc] 
1.0 kpc] 



(39) 



Using Zq — 0.1 kpc is an extreme assumption, which is 
certainly not realistic. Nevertheless, assuming values in the 
range 1-3 fxG seems reasonable. 
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Therefore, although there are uncertainties in the local 
value of the magnetic field, we suppose B ~ 1 /.jG in the fol- 
lowing, so that the ISRF model Ml complemented with the 
corresponding synchrotron losses leads to t; = 7.5 x 10^^ s 
in the Thomson approximation. The overall energy-loss rate 
b{E) /E"^ is plotted in the middle panel of Fig. [21 where it 
is shown to differ from the Thomson approximation with 
Ti — est — 10^^ s very often used in the literature, and 
which appears as the dashed straight line. In particular, 
we observe a cascade transition due to Klein-Nishina ef- 
fects, where it appears that the loss rate index a defined in 
Eq. Q decreases step by step from 2, its Thomson value: 
at about 1 GeV, relativistic corrections become sizable for 
interactions with the main UV component which is felt less 
and less by electrons; then, across the range 10-100 GeV the 
IR component gradually loses its braking potential, and fi- 
nally, above 10 TeV, interactions with CMB also cease. The 
value of the magnetic field sets the minimal value of the 
energy loss rate at higher energies. Since this latter is pro- 
portional to _B^, varying B from 1 to 3 translates into 
^ 1 additional order of magnitude in the energy-loss rate at 
high energy, as also depicted in the middle panel of Fig. [2] 
for completeness, we also display the case of taking 6 ^G. 
Note that considering CMB only provides a robust esti- 
mate of the minimal energy-loss rate, which converts into a 
maximal fiux by virtue of Eq. (j20p ; adding the synchrotron 
losses would instead define a next-to- minimal model for the 
energy losses. 

In the right panel of Fig. [51 we quantify the impact of 
using different energy-loss models to derive IS flux predic- 
tions, for which we adopt the med propagation setup and 
a template injection spectrum cx homogeneously dis- 
tributed in a thin disk. The dotted curve corresponds to 
the Thomson approximation with t; = lO^^s, where we re- 
cover a fiux with predicted index 7 = 7 + (5 + l)/2 = 2.85, 
as predicted from Eq. (|2ip with a — 2. The higher curve 
is the flux obtained with the minimal case for the energy- 
loss rate (the minimal t;), i.e. considering the CMB only, 
which provides the maximal flux. Indeed, in the Thomson 
approximation, the flux scales like ~ ^/rl =1/ -s/So i as seen 
from Eq. (|20p . The index reaches a plateau around 7 w 2.85 
in the range 10-1000 GeV, and then substantially hardens 
above 1 TeV because of relativistic effects. The next-to- 
minimal case exhibits the same feature, though the am- 
plitude is slightly reduced, as expected. Finally, we report 
the flux associated with our complete models Ml and M2, 
both associated with a magnetic held of 1 ^G (short dashed 
curves), 3 jiG (solid curves), and 6 ^G (long dashed curves). 
We remark that the naive prediction of 7 in the Thomson 
regime does not hold anymore, since the energy dependence 
of the energy loss a is no longer equal to 2, and the observed 
spectral index is signiflcantly harder. Indeed, we have to 
consider instead an effective value aes{E) < 2 to account 
for relativistic effects. Taking a larger value of the magnetic 
held slightly softens the index and decreases the amplitude, 
as expected. 

From this analysis of the local energy losses, we can es- 
timate that the related uncertainties translate into a factor 
of < 2 in terms of IS flux amplitude (</) cx y^), and, above 
10 GeV, ±0.1 in terms of spectral index (see the left panel 
of Fig. [21). Note, however, that this crude spectral analysis 
is valid for a smooth distribution of sources only, i.e. for 
secondaries. We see in Sect. 14.31 that considering discrete 



nearby sources of primaries strongly modifies this simplis- 
tic view. 

3. Secondary CR electrons and positrons 

We performed an exhaustive study of the secondary 
positron fiux in iDelahave et al.l ()2009() . which is qualita- 
tively fully valid for electrons and to which we refer the 
reader for more details. 

Secondary electrons originate from the spallation of 
hadronic cosmic ray species (mainly protons and a par- 
ticles) in the interstellar material (hydrogen and helium). 
This process produces also positrons, though different in- 
clusive cross-sections come into play. Since spallation in- 
volves positively charged particles, charge conservation im- 
plies that it gener ates more positrons than electrons {e.g. 
iKamae et al.ll20d6[) . This statement is not entirely accurate 
for neutron decay, but electrons arising from neutron decay 
have a very low energy (mostly E < 10 MeV), thereby out 
of the energy range considered in this paper. In sum, the 
steady-state source term for secondaries may in all cases be 
written as 

Q,{E,x)^A,:Y. J dE'ME',x)^^^^^^l^n,ix) , (40) 

where i flags the CR species of flux and j the ISM gas 
species of density n, the latter being concentrated within 
the thin Galactic disk, and d(Tij{E' , E) is the inclusive cross 
section for a CR-atom interaction to produce an electron 
or positron of energy E. 

For our default computation, we selected the proton- 
proton cros s- sectio n parameterizations provided in 
iKamae et al.l (|2006[ ). Any nucleus- nucleus cross-section 
{e.g. p — He or He — He) can be derived from the latter 
by applying an empirical rescaling, usually by means of 
a combination of the involved atomic numbers. However, 
this rescaling is found to be different for the production 
of 7r~ and tt"*", or e quivalently of e~ and e" *". We used the 
prescriptions from iNorburv fc Townsendl (|2007[ ) for this 
empirical rescaling. 

Fits of the proto n and a particle fluxes are provided in 
IShikaze et all (|2007[ ). based on various measurements at the 
Earth. Finally, we employed a constant density for the ISM 
gas, with uh = 0.9 cm~'^ and riHe = 0.1 cm~'^, confining 
these species to a thin disk of half-height h = 100 pc. This 
is summarized in cylindrical coordinates by 

n-iix) ^ e{h-\z\)e{R- r)nj . (41) 

For this spatial distribution of the gas, the spatial inte- 
gral of Eq. (ITOl) can be calculated analytically, following 
: Dclahave et a l. (2009); the solution is reported in Sect. lA.ll 
We underline that this approximation is locally rather good 
over the whole energy range as long as the true gas distri- 
bution does not exhibit too strong spatial gradients over 
a distance set by the half-thickness L — this is discussed 
in more detail for primaries in Sect. 14.21 In any case, this 
estimate is more reliable at high energies (^100 GeV) for 
which the signal is of local origin independent of L. At lower 
energies this approximation is expected to be valid for mod- 
erate L < A kpc, but much less trustworthy for large-halo 
models, as in the max propagation setup. For these extreme 
configurations, a more suitable description of the gas dis- 
tribution would be necessary. 
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Fig. 3. Flux predictions of secondary electrons (left) and 
positrons (right) at the Earth, for the min, med and max prop- 
agation models. 



In Fig. |3l we plot our results for the secondary electron 
(left-hand side) and positron (right-hand side) fluxes at the 
Earth. For the solar modulation, we used the fo rce-field a 
proximation with a Fisk potential of 600 MV (|Fisklll97ll 
The solid curves are derived with the Ml ISRF model, a 
relativistic tr eatment of energy l osses and nuclear cross- 
sections from iKamae et all (|2006D . The yellow band is the 
flux range available for all sets of propagati on parameters 
compa tible with B/C constraints derived in iMaurin et al.l 
(|200lt) . We observe that the min [max) configuration pro- 
vides the highest (lowest) and softest (hardest) flux, among 
the three beacon models. This can be clearly understood 
from Eq. (j20|) . since the min configuration is character- 
ized by the weakest value of diffusion coefficient normal- 
ization Kq associated with the strongest index 5, which is 
reversed in the max configuration. However, we recall that 
these models were named so for sources distributed all over 
the diffusion halo, not only confined to the disk as is the 
case here. 

The discussion in IDelahave et al] (120091) on the theoreti- 
cal uncertainties affecting secondary positrons is fully valid 
for secondary electrons. Aside from energy losses, errors 
may originate from either uncertainties in the light nuclei 
flux, or uncertainties in nuclear cross-sections, or both. The 
former can be evaluated by using different fits of the local 
measurements, and by retro-propagating the CR nuclei fiux 
to account for potential spatial gradients. The latter may 
be estimated by considering alternative parameterizations 
of nuclear cross-sections. This is illustrated in Fig. [3] with 
the short- dashed curves co mputed with the nuclear cross 
sections of lTan fc Nd (I1983D , which are shown to differ from 
our default model only at low energy below a few GeV. All 
these effects were studied for positrons in IDelahave et al.l 
(j2009t) and lead to an uncertainty of about 40 %; this theo- 
retical error is also valid for secondary electrons. Note, how- 
ever, as we see later, that the electron flux is most likely to 
be dominated by the primaries, in contrast to positrons for 



which secondaries are a major component. Therefore, un- 
certainties in the secondary contribution has more impact 
for positrons than for electrons. 

Finally, we emphasize that the present re s ults d iffer 
slightly from those derived in IDelahave et al.l (|2009l ) be- 
cause the energy losses are now treated in a fully relativistic 
formalism. Not only does this slightly change the normal- 
ization at low energy by a factor ti /tdo9 ~ 0.9, but, more 
importantly, this hardens the spectral shape due to Klein- 
Nishina effects. This is more striking in Fig. [31 where the 
long-dashed curves are the predictions calculated in the med 
configuration and the Thomson limit for the energy losses. 
In Appendix [B] we provide user-friendly fitting formulae 
that closely reproduce our calculations of the secondary 
electron and positron fiuxes in the med propagation setup. 



4. Primary electrons 

The GeV-TeV CR electron flux at the Earth is dom- 
inated by a primary component originating from elec- 
trons accele rated in both SNRs and pulsar wind neb- 
ula3(PWN) (|Blandford fc Eichled [T981 . CR sources are 
therefore connected to the explosion of supernovae(SNe), 
and we discuss a more general framework in Sect. [5l 
Predicting this primary contribution is a rather difficult ex- 
ercise because it involves characterizing the energy distribu- 
tion of these electrons at sources and their spatial distribu- 
tion, in addition to their transport to the Earth. Moreover, 
since GeV-TeV electrons have a short range propagation 
scale, local sources are expected to play an important role. 
Rephrased in statistical terms, since the number of sources 
exhibits large fluctuations across short distances, the vari- 
ance affecting the predictions is expected to increase with 
energy, when the effective propagation volume decreases. 
Therefore, fluctuations in the properties of local sources 
will have a strong impact. In contrast, the low energy part 
of the CR electron spectrum might be safely described in 
terms of the average source properties, namely a smooth 
spatial distribution associated with a mean injected energy 
distribution. 

In the following, we estimate the primary flux of elec- 
trons and quantify the associated theoretical errors. More 
precisely, we wish to quantify the relative origin and im- 
pact of these uncertainties. To do so, we first discuss in 
Sect. 14.11 the spectral shape properties that we consider in 
the forthcoming calculations, focusing on SNRs for the mo- 
ment (pulsars are discussed in Sect. |5l together with the 
primary positrons); available spatial distributions are pre- 
sented in Sect. 14.21 Then we discuss the uncertainties asso- 
ciated with the modeling of a single source in Sect. 14.31 We 
finally discuss the primary fiux and related uncertainties in 
Sect. 14.41 in which we make a thorough census of the local 
SNRs likely affecting the high energy part of the spectrum. 



4.1. Spectral properties of SNRs and related constraints 

Most SNR models (e.g. lEUison et al.ll2007l: [Tatischefai2009t) 
rely on the acceleration of CRs at non-relativistic shocks 
fe.o. lMalkov fc O'C Drurvll200lh and predict similar energy 
distributions for the electrons released in the ISM, which 
can be summarized as 



E 



(42) 
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The spectral index 7 is usual ly found to be ar ound 2 
over a significant energy range (jEllison et al.l[2007t ) — but 
to exhibit large variations at the edges — in agreement 
with radio observations. Gamma-ray observations suggest 
that the energy cut-off Ec is greater than a few TeV {e.g. 
lAh aronian et al.|[2009bf ). These studies find rather similar 
indices for protons and electrons. Since protons of energy 
above a few GeV are barely affected by energy losses and 
have a long-range propagation scale, the proton spectrum 
measured at the Earth can provide information about the 
mean index at sources. Since 0p oc Q{E)/ K{E) oc E^"' , the 
index at source is therefore 7 « 7 — (5. With 7 ~ 2.8 and 
5 in the range 0.5-0.7, one finds 7 in the range 2.0-2.3, in 
rough agreement with theoretical predictions. We discuss 
complementary constraints from radio observations at the 
end of this subpart. 

Aside from the spectral shape, sizing the value of the 
normalization Qo is much more problematic. To describe a 
distribution of sources in the Galaxy, one usually assumes 
that the high energy electron injection is related to the 
explosion rate of SNe, so that we can suppose that Qo is 
such that the total energy carried by electrons is given by 



dE'E' Q{E')^fE^T^ 



(43) 



where F^, is the SN explosion rate, E^, is the kinetic energy 
released by the explosion, and / is the fraction of this energy 
conferred to electrons. Since we are only interested in the 
non-thermal electrons, we assume that E^min = 0.1 GeV. 
Note that the spectral index influences the normalization 
procedure sketched above. For a single source, the same 
expression holds but without F* — we discuss this case 
later on. For convenience, we further define the quantities 



= fE, , 
F_B = T-i^E-), ^ 



(44) 



which also helps us to discuss the normalization issue. 

Constraining F^,, E^,, and /, is a difhcult exercise. The 
explosion rate of SNe is typic ally predicted to be a 1-5 per 
century and per galaxy {e.g. Ivan den Bergh fc TammannI 
Il99lt iMadau et al.l 19981). which is consistent with obser - 
vations (e.g. IValinia fc Marshalll 119981: iDiehl et al.l |2006[) . 
Nevertheless, SNe are of different types, and may thereby 
lead to different CR acceleration processes. About 2/3 of 
SNe are expected to be core-collapse SNe (CCSNe), the 
remaining 1/3 consisting of type la SNe (SNela). This is 
at variance with the statistics derived from observations, 
which identify a higher percentage of the more luminous 
latter type. CCSNe are often observed in star formation re- 
gions and spring from the collapse of massive stars > 8 Mq , 
while SNela, produced by older accreting white dwarfs, are 
more modest systems, while being more widely distributed. 

Explosions of CCSNe with masses < 2QMq can typ- 
ically liberate a huge amount of energy, ~ iQ53-5i ^^^^ 
~ 99% of which is released in the form of neutrinos 



neutron star {e.g . iBurrow I 12OOOI: IWooslev fc .Tml 



e proto- 
fe] 120051: 



when integrated after the cooling phase of the proto- 
a r (e.g . 

iJanka et aLll2007l) . The total kinetic energy available from 
CCSN explosions, mostly due to the energy deposit of the 
neutrinos in the surrounding material, is about ~ 10^^ 
erg. CCSNe usually give rise quite complex systems char- 
acterized by SNRs beside (or inside) which one can find 



active neutron stars such as pulsars and associated wind 
nebula (PWN) — we focus on pulsars in Sect. [5j In con- 
trast, SNela are much more modest systems, the explo- 
sions of which blast all the material away without giv- 
ing birth to any compact object and release an energy of 
about 10^^ erg, which is transferred entir e ly to the sur- 



rounding me dium {e.g. Nomoto et al.l 119841 : iGamezo et al.l 



i2003.: Mazzali et al.|[2007 ). It is interesting to remark that 
despite the quite different natures of CCSNe and SNela, 
both types feed the ISM with the same typical kinetic en- 
ergy. Since we consider only SNRs in this part, we assume 
that Ei, = 10^^ erg in the following. 

The frac tion of SN energ y conferred to electrons was 
studied in lTatische"3 ()2009r ). and found to be about 
/ ^ 10~^ — 10~^. This result is rather independent of 
the exact values of the spectral index 7 and the cut-off en- 
ergy Ec, but the theoretical error is still of about one order 
of magnitude. 

At this stage, we emphasize that the theoretical uncer- 
tainty in Te = f E^,^^ already reaches about 2-3 orders of 
magnitude on average, which is quite huge and translates 
linearly in terms of flux. 

For known single sources, it is possible to derive tighter 
constraints on the individual normalizations Qo from ob- 
servations at wavelengths for which electrons are the main 
emitters. This is precisely the case for the non-thermal ra- 
dio emission produced by synchrotron processes, provided 
the magnetic field is constrained independently. The syn- 
chrotron emissivity associated wit h an electron source of 
injection rate Q{E) is given by (iGinzburg &: Svrovatskiil 



119651: iBlumenthal fc Gouldll 19701: lLongaidll994D 



J{i^) = \J^dd sm{9) J dvPs{v,e) Q{E = hv) , 



(45) 



where an average is performed over the pitch angle 0, and 
where the synchrotron radiation power is defined as 



Ps{i^,e) 

where x 



X / dyK^i^{y) 

V 3 2 3 6^7^ 

— ; vs = -i — , 

J/, 2 4 TT m 



(46) 



and Va, which is called the synchrotron peak frequency, cor- 
responds to the average frequency of the synchrotron emis- 
sion arising when an electron of Lorentz factor 7 interacts 
with a magnetic field B, and is cx 7^t'c, where Vc is the 
cyclotron frequency. We note the parallel between the syn- 
chrotron process and the the inverse Compton process: for 
the latter, an initial photon of energy Ek would indeed be 
boosted to ^ 7^i?/c- We could calculate the radio flux from 
the emissivity, but it is also interesting to derive a more 
intuitive expression, wh ich actually provides a fair approx- 
imation (|Longaiijll994f ) 



d(l3{l^) 

dv 



dv : 



Q{E) dE , 



(47) 



where we assume that the entire radiation is emitted at the 
synchrotron peak frequency v = Vs, which links v to E. The 
distance from the observer to the source is denoted by d. 
We remark that only the synchrotron part of the electron 
energy loss rate b{E) appears — we therefore assume that 
this is the most efficient process within the source — and 
that we neglect the possible reabsorption of the synchrotron 
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emission. Since b{E) is the energy lost by an electron, it 
corresponds to the energy of the emitted photon, so the 
factor l/{hv) allows us to infer the number of photons. 

This expression allows us to constrain Qo by means of 
the source radio brightness Briiy), which is usually found 
in catalogs 



1 



dv' 



(48) 



sync ^/pN dE 



We readily derive 



(49) 



which translates into 

1.2 X 10''^ X (0.79)'^ 



GeV- 



(50) 



kpc 



GHz 



B 



100 AtG 



Jy 



We have just recovered the well-known relation between the 
radio index and the electron index, 7r = (7 — l)/2. 

An up-to-date catalog of SNRs can be found in iGreenI 
(|2009D . which contains ~265 objects. Among these objects, 
only 70 have estimated distances to the Earth, and 207 have 
measured radio spectral indices. Observations, however, are 
not expected to reflect the actual statistical properties of 
the whole population of Galactic SNRs because of obser- 
vational selection effects favoring the brightest sources and 
sites of fainter background (high longitudes, towards the 
anticenter). Disregarding the spatial distribution of these 
objects, which is probably strongly biased, this sample may 
still b e fairly repre sentative of their general spectral prop- 
erties (|Greenll2005l) . 

We compiled histograms of the measured radio in- 
dices and the estimated intrinsic luminosities — £/(4 tt) = 
d^i3i.(lGHz) — in the right and left panels, respectively, 
of Fig. m The radio indices clearly appear to exhibit a 
Gaussian distribution, whereas luminosities follow a log- 
normal distribution. This points towards similar physical 
grounds for the electron properties at sources, which is 
obviously unsurprising. With these distributions, we can 
derive mean values and statistical ranges for the parame- 
ters. We find that (7^ = 0.50 ±0.15 and {d^ Br{l GHz)) = 
exp {6.26 ± 1.95} Jy.kpc^. One can therefore infer that the 
electron index (7) = 2 (7,.) -|- 1 = 2.0 ± 0.3 in very good 
agreement with theoretical expectations. Although this re- 
lation between the radio index and the electron index is not 
entirely accurate (because of other radio components or ab- 
sorption), and although some systematic errors also affect 
the data, this provides a complementary means of sizing 
the uncertainty, which is consistent with that of theoretical 
results. 

We use this statistical information to directly constrain 
the single source normalization Qq from Eq. (|50p . but we 
need to estimate the magnetic field in SNRs. From the ob- 
servational point of view, information about the electron 
density and magnetic field at sources is degenerate. More 
insights may come from theoretical studies of the amplifi- 
cation of magnetic fields in sources from numerical simula- 
tions, which involve CRs themselves as seeds and amplifiers. 



iLucek fc Bell 
with observa- 
we finally find 
which trans- 



lO''^ erg (with 
^reement with 
biased, as ex- 



The current state-of-the-art simulations {e.g. 
I2OOOO support B - 100 AiG, in agreement 
tions, and that we use here. With this value, 
(Qo) = 3.9 X lO''^ GeV-i for an index 7 2 
lates into (E^) ~ 4.3 x 10^° GeV ~ 6.9 x 
a cut-off Ec ~ 10 TcV). This is in rough a^ 
the other values derived above, but probably 
pected, towards the brightest objects. 

4.2. Spatial distribution of sources 

Although GeV-TeV electrons have a short-range propaga- 
tion scale, the injection rate of energy discussed above is 
insufficient to describe the Galactic CR electrons. We need 
to specify the spatial distribution of sources. For nearby 
sources, for which observational biased are less prominent, 
we can use available catalogs, which may provide a rather 
good description of the local CR injection. Nevertheless, 
for more distant sources, which have influence on the in- 
termediate energy range ~1-100 GeV, we have to rely on a 
distribution model. 

Since 2/3 of SNc are expected to be CCSNe, one can 
use pulsars as tracers of the SNR distribution, instead 
of SNRs themselves, the observed population of which is 
much more modest. As an illustration, the ATNF cata- 
log (jManchester et all 120051 ) lists more tha n 1800 pulsar s 
compared to the 265 SNRs contained in iGreenI (|2009[ ). 
Nevertheless, a too naive use of the statistics would lead to 
errors since it is well known that data do not reflect reali ty 
faithfully because of detection biases (e.g. iLorime H I2004D . 

There are few distribution models available in the lit- 
erature. Since the energetics associated with the source 
injection (birth) rate has been discussed above, we are 
only interested in the normalized source distribution here. 
Consequently, the normalization coefficient in front of each 
model is fixed such that it normalizes the spatial distribu- 
tion to unity within the diffusion halo characterized by its 
radius R and half-thickness L. Moreover, in the following, 
we set the position of the Sun at Rq = 8 kpc from the 
Galactic centeiQ. 

Most of models have radial and vertical dependences of 
the form 



p{r, z) = pq r° exp 



ro 



exp 



(51) 



where po ensures the normalization to unity. For simplic- 
ity, we discuss only differences of the radial distributions in 
the following, since the vertical distribution is fairly sim- 
ilar among studies. We thereforee keep fixed the vertical 
dependence as in the above equation, with zq = h ~ 0.1 
kpc, throughout the paper. 

Different sets of values c a n be found in the literature 
for the pair (a, rn). Lorimer (f2004D her eafter L04, found 
(2.35, 1.528 kpc): lYusifov fc KiicaH(l200l . h ereafter YK04, 
derived (4, 1.25 kpc): while lPaczvnskil ()1990D . hereafter P90, 
early determined (1, 4.5 kpc). Finally, in contrast to the pa- 
rameterizati on sketched above, we recall the distribution 
proposed bv lCase fc Bhattacharval (|1998f ). hereafter CB98, 



^ |http : / / www ■ at nf . csiro . au/research/pulsar/psrcat 
^ Some of the distributions listed in this paragraph arc ac- 
tually derived assuming 8.5 kpc, but we disregard this small 
change to make the discussion easier. 
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though it was obtained from a fit to data of poor statistics 
for 36 SNRs 

p(r,z) = po sm(^TT^+e^ cxp|-^| exp|-M| ^ (52) 

where we have added the same vertical term as in Eq. (|5ip . 
The authors found ro = 7.7 ± 4.7 kpc, = 17.2 ± 1.9 
kpc, and 6 = 0.08 ± 0.33. This relation is only valid for 
r < fsfl — ^/tt ), i.e. within 1 6 .8 kpc , and null beyond. Note, 
however, that iBrogan et all (I2006D reported the detection 
of 35 new remnants in the inner Galaxy, and suggest that 
former radial distribution estimations should be revised. 

To understand the deviations induced in the electron 
flux prediction when using these different distributions, it 
is convenient to define the following halo function 

Tq(A) = J d^Xs p{rs, Zs) g\{X, Xq ^ Xs) (53) 
where C/a(A, ajg a;,,) = b{E)g{E,XQ ^ Es,Xs) , 

which determines the probability of an electron reaching 
the Earth given its propagation scale A — see Eq. ^ — 
and the normalized spatial distribution of source p. The 
electron flux is the the energy integral of the product of this 
probability and the source spectrum, such that the shape 
of this probability function provides a preliminary taste of 
the final result. More importantly, it allows us to connect 
the spatial origin of the signal with energy, through the 
propagation scale A. 

In the left panel of Fig. [SJ we plot the spatial distribu- 
tions listed above as functions of the galactoccntric radius 
r, and in the galactic plane (z = 0). We see that except for 
the solar neighborhood, where relative amplitudes can vary 
by a factor of ^2 at most, the spatial distributions in the 
direction of the Galactic center and towards the anticenter 
are quite different from each other. Nevertheless, these dif- 
ferences are significantly lower in terms of Tq(A), because 



of the spatial average — see Eq. ([55)1 . This is shown in the 
middle panel of Fig. [SI where we have plotted A x Tq (A) as 
a function of A for the different spatial distributions and for 
the mm, med, and max propagation setups. We see that the 
probability is maximal and constant — A Tq grows linearly 
with A — for short propagation scales up to X h = 0.1 
kpc. Then, the probability decreases linearly with A — AT 
exhibits a plateau — before shrinking exponentially when 
X ^ L, L being larger and larger from the min setup to the 
max setup. Each spatial distribution model is characterized 
by a very similar curve that differs mostly in terms of am- 
plitude. This can be understood in the following manner: 
when X < h, the source can be considered as homogeneous 
in 3D space, then Tq oc / drr'^X~^ exp{— r^/A^} = est; 
when h < X < L, since the source distributions do not 
exhibit strong radial variations on the kpc scale, they can 
be considered as thin disks, and one recovers the solution 
T cx A~^ derived in Eq. ((20l) : for A > L, electrons escape 
the diffusion zone. This points towards the possibility of 
modeling, while only locally, the source distribution with 
a z-exponential infinite disk, for which full analytical solu- 
tions of the spatial integral exist. The green curve in the 
middle panel of Fig. [5] is the z-exponential disk approxima- 
tion associated with P90, as an illustration, and is shown 
to provide a rather good approximation except for large 
diffusion thickness L > 4 kpc. The z-exponential disk ap- 
proximation is defined in cylindrical coordinates as 

Pd(r,z)=p(i?0,O)0(i?-r)exp|-^| , (54) 

where p(i?0,O) is the local value of the normalized den- 
sity given in Eq. ([5T|) . This approximation is valid for lo- 
cal predictions provided the spatial distribution p does not 
vary significantly over a distance ^ L, which is the case 
for moderate L. In the right-hand side panel of Fig. [SJ we 
compare the disk approximation with the full calculation in 
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Fig. 5. Left: Spatial distribution models for SNRs and pulsars. Middle: Corresponding halo functions defined in Eq. (|53|) and that 
characterize the transport probability as a function of the propagation scale A, which decreases with energy. Right: Ratio of flux 
predictions to their associated exponential disk approximations. 
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Fig. 6. Fraction of the signal reaching the Earth as a function 
of the integrated radius, for different energies and different spa- 
tial distribution models, using full relativistic energy losses. The 
Thomson approximation result is reported in dashed line. 



terms of fluxes: for difFcrcnt spatial distributions, we plot 
the ratio approximated flux / exact flux for our three bea- 
con propagation setups. We can sec that the exponential 
disk approximation is quite good above a few GeV for the 
min and med cases, as expected, having an accuracy better 
than 5%. Errors are obviously larger in the max case be- 
cause of the larger spatial gradients exhibited by the spatial 
distributions within L ~ W kpc. 

A final useful exercise regarding the smooth spatial dis- 
tribution modeling consists of checking the cumulative frac- 
tion of the IS signal received at the Earth as a function of 
the radial integration distance. In Fig.|6l we report this frac- 
tion for spatial model L04 at different energies, assuming an 
injection spectrum cx E~'^, and for both the Thomson ap- 



proximation and the relativistic energy losses. We see that 
this fraction increases more quickly at high energy than at 
low energy, as expected from en ergy losses. This is con sis- 
tent with the result obtained in iDelahave et al.l (|2009t) for 
secondary positrons. Nevertheless, above '--^lO GeV, we can 
observe that relativistic effects come into play and a dif- 
ference appears between the Thomson approximation case 
and the relativistic case. Indeed, the latter induces a longer 
propagation scale at high energy, and consequently softens 
the rise of the cumulative fraction. This would be slightly 
less significant for a magnetic field of 3 instead of 1 fiG, 
though still observable. 

Another important piece of information that we can 
derive from Fig. [5] is that the cumulative signal fraction 
is > 95% (80%) for r > 2 kpc and E > 100 (10) GeV. 
This helps us to define consistent means of including local 
sources in our predictions, as we discuss later in Sect. 14.4.21 
Indeed, we know at present that if we replace the smooth 
spatial distribution within 2 kpc with discrete sources, these 
latter can affect the whole available energy range quite sig- 
nificantly: if powerful enough, local sources will dominate 
above a few tens of GeV, otherwise, flux predictions will be 
significantly depicted compared to a smooth-only descrip- 
tion of sources, for a given normalization pattern. 

4.3. Sizing the uncertainties for local sources 

Before discussing the contribution of local known SNRs to 
the GR electron flux (see Sect. I4.4.2p . it is essential to re- 
view the impact of uncertainties in the main parameters 
describing the source. They are only a few, but their effects 
on the flux are shown to be important and degenerate. 

Apart from the propagation modeling and related pa- 
rameters that were presented in Sects. 12.21 and 12.51 theo- 
retical errors may originate from uncertainties (i) in the 
spectral shape and normalization, (ii) in the distance esti- 
mate, (iii) in the age estimate and (iv) in our understanding 
of the escape of cosmic rays from sources. The last point is 
actually still debate d and poorly known in detail (see e.g. 
iGaprioli et al.l[2b09t) . though it is clear that the release of 
cosmic rays in the ISM is a time- and energy-dependent pro- 
cess which takes place over ~ 10^~^ yr, i.e. the lifetime of 
the source. Since this timescale is still almost always much 
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Fig. 7. Main uncertainties associated with the flux of primary electrons injected from a single bursting source. Top left: varying 
the age at a fixed distance of 500 pc — notice that taking 1 kpc would have suppressed the 3 kyr source for causality reasons. Top 
right; varying the distance for a fixed age of 5 kyr. Bottom left; varying the spectral index for fixed age of 5 kyr and fixed distance 
of 500 pc. Bottom right; varying propagation parameters, with the min, med and max setups from Table [1] 



lower than the diffusion timcscale, — cP/^Kq ~ 0.1 — 10 
Myr for distances in the range 0.1-1 kpc, we ignore the dy- 
namical aspects of injection in this study, while we stress 
that they may lead to sizeable effects, especially in the case 
of very nearby sources. The first point was discussed in 
Sect. 14. Il and is featured by two main parameters: the spec- 
tral index at source 7 and the energy released in the form of 
high energy electron / E^,, both related in the normalization 
procedure given Eq. that allows to derive Q/. Points 
(ii) and (iii) have some impacts that can be understood 
from Sect. 12.21 and Sect. 12.31 Although the consequences of 
varying these parameters can be understood from equations 
only, we aim here to illustrate them in a more pedagogical 



way. To do so, wc will consider a template event-like source 
located in the Galactic plane (2 = 0) at a distance d to the 
Earth and bursting a population of electrons a time (age) 
ago: 

Q^{t„Es, X,) = S{t, - U) 5{z,) 5^{r, - d) Q{E,) , (55) 

where the spectrum Q{Es) is given by Eq. (|^^ . Wc will 
assume here that / iJ^ = 2 x 10''^'' erg. Note that, as empha- 
sized in point (iv) above, a more realistic source term would 
not involve a burst-like release of electrons in the ISM at 
time t^, but instead a more complex time-dependent en- 
ergy spectrum. Such refinements arc beyond the scope of 
this paper. 
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In Fig. [71 we plot the electron flux for different config- 
urations of the parameters, the default configuration being 
defined by: med, 7 = 2, E'^ = 10 TeV. 

In the top left panel (a), we show the source age ef- 
fect; in the top right panel (b), we illustrate the distance 
effect; in the bottom left panel (c), we sketch the spectral 
index effect; while in the bottom right panel (d), we plot 
the propagation model effect. For all panels, we report the 
fluxes calculated in both the Thomson approximation and 
the full relativistic treatment of the energy losses, as dis- 
cussed in Sects. and 12.51 

As a first comment, we emphasize that the Thomson ap- 
proximation can lead to a very strong under-estimate of the 
spectral break inferred from energy losses, up to one order 
of magnitude in the examples shown. This is a mere conse- 
quence of the over-estimate of the energy loss rate at high 
energy. The net effect obviously depends on the magnetic 
field and on the actual cut-off considered at the source. As 
regards the latter, we see that using a value of 10 TeV al- 
ready induces an underestimate by a factor of ^5-10 of the 
break predicted in the non-relativistic regime. Many studies 
of the topic have employed the Thomson approximation. 

The second important comment to make is that it is ac- 
tually quite difficult to relate the observed spectral index to 
the source spectral index, because of the complex and de- 
generate effects coming from all parameters: distance, age, 
source index, energy cut-off, normalization and diffusion 
coefficient. For instance, we see that a large diffusion coef- 
ficient (mm model) can make a source of 300 kyr resemble 
a source of 30 kyr associated with a larger diffusion coeffi- 
cient and a lower energy cut-off. In any case, a mere glance 
at the four panels of Fig. [7] is striking enough. 

This exhaustive analysis of the impact of the main pa- 
rameters characterizing individual sources already points 
towards the difficulties that we encounter in the interpre- 
tation of the data. Nonetheless, although this part might 
look depressing at first sight in the perspective of making 
predictions, it is still very useful to estimate the theoretical 
confidence level of our forthcoming attempts. 

4.4. Primary electron flux and theoretical uncertainties 

In the previous parts of this section, we have discussed the 
main physical quantities relevant for predictions of the pri- 
mary electron flux at the Earth, emphasizing their role as 
potential sources of uncertainties. Here, we implement the 
full calculation and compare our results with available data 
on the electron flux. We stress that pure electron data are 
not numerous and rather old, since most of recent exper- 
iments either do not distinguish electrons from positrons 
or have not yet released their charge-discriminating data. 
We therefore only use the electron d ata from CAP RICE 
(fBoe zio e t al]|2000D, HEAT (iPuVcrn ois et all |2001[ ). and 
AMS-01 |Alcaraz et al.|[2000l ). to avoid any confused inter- 
pretation mixing positrons. The pure positron case and the 
full case are discussed in Sects. [S] and [51 respectively. 

We first compare the predictions arising from a smooth 
description of sources, for which we adopt the L04 spa- 
tial distribution. Indeed, we demonstrated in Sect. 14. 2 1 that 
using different spatial distributions causes only small dif- 
ferences in the overall flux normalization locally. 

Wc then estimate the contributions of all known local 
SNRs that can be added to a smooth and more distant com- 



ponent, following the method proposed in lKobavashi et al.l 
(l200l . 



4.4.1. Smooth description of sources 

Our model for a smooth distribution of SNRs includes a 
propagation setup, a spatial distribution (here L04) and an 
injected spectrum, and it is interesting to check some of 
the possible configurations against the data. In particular, 
we attempt to constrain the injection normalization neces- 
sary for a model to fit, at least roughly, the data. For the 
spectrum, we test different spectral indices, but keep the 
energy cut-off at 3 TeV. As a reference normalization, we 
use a SN explosion rate of 4/century, a SNR total energy of 
E^, = 10^^ erg, of which a fraction of / = 2 x 10~^ is carried 
by electrons, giving therefore F^ = 8 x lO''^ erg/century. 

In Fig. [51 we report various flux calculations, for which 
we applied a solar modulation correction with a Fisk po- 
tential of 600 MV. In the left panel, we show the effect of 
varying the injected spectral index from 2 to 2.4 for the 
med propagation setup, using both the Thomson approxi- 
mation and the relativistic regime for the energy losses. In 
this plot, we have renormalized F^; by a factor of 5 for all 
indices, so that we see that reasonable fits to the data can 
be obtained within the expected normalization range dis- 
cussed in Sect. 14.11 This means that the expected energy 
budget available for electrons is in rough agreement with 
what is needed to explain the current observations. From 
the same plot, we could also conclude that the injection 
spectral index should be slightly softer than 2. Nevertheless, 
this also depends on the logarithmic slope 6 of the diffu- 
sion coefficient, as seen from Eq. (j2ip — complementary 
constraints on j + S could also be derived from high en- 
ergy proton data, based on the assumption that the proton 
index is the same as the electron index after their accel- 
eration at sources and that proton propagation is simply 
described by diffusion {i.e. neither reacceleration nor con- 
vection) . This is illustrated in the left panel of Fig. [51 where 
we show the effect of the theoretical uncertainties in the 
propagation parameters, using the same spectrum normal- 
ization and the same spectral index for all models. We see 
that the min model gives the larger amplitude because of 
its smaller value of Kq and the softer observed index due 
to its larger diffusion slope 6 (see Table [T] and Eq. (I^Uj) ) 
— the analysis is reversed in the max configuration. For a 
given normalization, the amplitude uncertainty is therefore 
proportional to \/Ko, which gives a factor of ~ 7 from the 
min to the max configurations. In both panels o f Fig. [51 we 
also re port the pre diction obtained inlMoskalenko fc Strong 
(I1998D . as fitted in lBaltz fc Edsiol (Il998l ). where the authors 
used an injection index of 2.1 below 10 GeV, steepening to 
2.4 above. This model, very often quoted as a reference 
model, is shown for comparison. 

It is noteworthy that since the data have a quite limited 
statistics and range up to ^ 40 GeV only, they are prob- 
ably insufficient to provide strong constraints on the elec- 
tron cosmic ray component. Moreover, we recall that this 
smooth description of the SNR contribution is not valid lo- 
cally above a few tens of GeV, where we expect discrete 
effects to become important. Nevertheless, this preliminary 
analysis is still useful to delineating the relevant ranges of 
the spectral index and the injected energy. Likewise, it helps 
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Fig. 8. primary electron flux from a smooth distribution of SNRs. Left: fluxes associated with injection indices 7 from 2.0 to 2.4. 
Right: propagation effect for 7 = 2.3. A renormalization factor of 5 has been applied to Fb in both panels and a cut-off energy of 
3 TeV has been considered. 



us to determine the influence of distant sources relative to 
local ones. 

Finally, we emphasize that we only considered a sin- 
gle contribution from a SNR population. Nevertheless, the 
electron-positron pair injection from pulsars is also likely 
to account for a significant additional contribution to the 
local electron budget. It is not clear whether this contri- 
bution should have the same spectral index, and one could 
for instance model the smooth electron component with 
a combination of two spectral components, leading to an 
additional freedom in the normalization procedure. This 
electron component from pulsars is discussed in Sect. [5j 

4.4.2. Contributions from known local sources 

As discussed earlier, contributions from local sources arc ex- 
pected to be significant a bove a few tens o f GeV. Following 
the method proposed in iKobavashi et al.l (|2004[) , we take 
a census of all known sources of primary electrons located 
within 2 kpc from the Earth in order to compute their as- 
sociated fiux explicitly. 

To proceed, we first took adva ntage of the informa- 
tion provided in the SNR catalog of iGreerJ (j2009t l. and we 
performed an extensive synthesis of all published proper- 
ties and associated errors (mostly from radio data). We 
found 26 SNRs within 2 kpc in this catalog, to which 
we added an extra-obje ct, Antlia SNR (jMcCuUough et all 
[200l lSh inn et alll2007t) . A full description of these sources 
including information about distance, age, spectral index, 
radio flux, associated objects, and bibliographic references 
is available in Appendix [Cj These properties are summa- 
rized in Table ICTl 

The event-like contribution of a single source is read- 
ily computed from the results obtained in Sect. 12.21 As a 
word of caution, however, we stress that the time argu- 
ment used to feed the time-dependent propagator given in 





Ec 


1 TeV 


10 TeV 


100 TcV 


e 
S 


7 = 2.0 


5, 20, 22, 24 


+ 23 


+8, 18, 19, 26 


7 = 2.2 


5, 20, 22, 24 


-f 8, 18 


+ 17,19 


7 = 2.4 


5, 20, 22, 24 


+8, 17, 18 


+ 19 


med 


7 = 2.0 


5, 20, 22, 24, 
26 


+ 4, 11, 19 


+8, 18, 23 


7 = 2.2 


5, 11, 20, 22, 
24 


+4, 8, 18,23 


+ 19 


7 = 2.4 


5, 11, 20, 22, 
24 


+4, 18, 19, 23 


+ 8, 17 


max 


7 = 2.0 


5, 20, 22, 24, 
26 


+ 8, 11, 19, 
23 


+4, 18 


7 = 2.2 


5, 20, 22, 24 


+8, 11, 18, 19 


+ 4, 23, 26 


7 = 2.4 


5, 20, 22, 24 


+8, 11, 18, 19 


+4, 23, 24 



Table 3. Ranked id numbers (see Table lC.lj) of those SNRs 
that contribute more than 10% of the signal for various 
propagation models, cut-off energies ec, and spectral index 
7. Index 7 is only used for sources that are not constrained 
by radio data. The dominant sources are the Cygnus Loop 
(5), Monogem Ring (20), Vela (22), and Antha (24). 



Eq. should not be the observed age of the object given 
in catalogs, but instead the actual age, equal, in principle, 
to the observed age plus d/c. Indeed, most of the age es- 
timates depend on the dynamical properties of the objects 
inferred from multiwavclength observations, which corre- 
spond to the properties the object had a time d/c ago. For 
the injection spectrum, we utilize Eq. and set the spec- 
tral index 7 from the observed radio index 7^ — 7 = 27^-!-!. 
We constrain the spectrum normalization with the observed 
radio fiux using Eq. ([50)) . 

Although SNRs are expected to provide an important 
contribution to the primary electron flux, we emphasize 
that pulsars are also expected to produce and acceler- 
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Fig. 9. Left: Plot of the observed age versus distance to the Earth for our sample of local SNRs (and associated uncertainties, 
see Table [CTTjl . The dashed lines correspond to limits beneath which a local source cannot contribute significantly to the signal at 
the corresponding energy (valid only in the med propagation model — see Table [l]). Indeed the age sets an upper limit, while the 
distance sets a lower limit to the energy range - see Sect. 14.41 Right: Same plot for our complete sample of local SNRs and pulsars. 
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Fig. 10. Primary electron flux from local SNRs in the med propagation model and using radio observation constraints. Left: 
energy cut-off at £c = 1 TeV. Middle: Er_ = 10 TeV. Right: = 100 TeV. 



ate electron-positron pairs. Modeling the electron injection 
from pulsar is discussed in more detail in Sect. [SJ to which 
we refer the reader. We found that ^ 200 pulsars located 
within a distance of 2 kpc from the Earth could contribute 
to the local electron budget, among which few may be dom- 
inant (see Sect. 15. 2p . For consistency reasons, we have to 
include the contribution of these pulsars to the electron 
flux. 

An additional important remark should be made about 
the non-observed local sources of primary electrons that 
should exist. So far, we have listed 27 SNRs and about 200 
pulsars. Nevertheless, we recall, as we discuss in more de- 
tail in Sect.[Sl that pulsars arc rotating neutron stars origi- 
nating from core-coUapsc supernova explosions. Therefore, 
each pulsar should be accompanied by a SNR. Such a sys- 
tematic association is obviously not supported by observa- 
tions. This is already illustrated in our object list, in which 
we find only 27 SNRs for 200 pulsars. Among these 27 lo- 
cal SNRs, only 10 have a known pulsar counterpart (very 
often with differences in their distance and age estimates). 



However, this certainly does not mean that the theoreti- 
cal expectation is wrong, since not all SN explosions lead 
to pulsars, but instead that the counterparts are proba- 
bly not bright enough to be observed. Therefore, again for 
theoretical consistency reasons, we choose here to add a 
SNR counterpart to each non-associated pulsar, but with 
a brightness such that it could not be observed with cur- 
rent telescopes. We adopted B{1 GHz) < IJy as a general 
criterion for non-observed SNRs. These local statistics can 
be tested against predictions of the SN explosion rate T^,. 
If we assume at zeroth order that sources are distributed 
homogeneously inside a flat disk of radius i? = 20 kpc, then 
the local explosion rate within a radius of ?■ = 2 kpc around 
the observer is given by rj,(r/i?)^ = 0.01 x F^,, leading to 
~0.01 SN/century for usual values of F^,. This can be com- 
pared with the observed local explosion rate, which we can 
estimate from the number of sources in our sample divided 
by the oldest age, i.e. ~ 200/30 Myr ~ 0.7 x IQ-^/century. 
This rough calculation leads to a difference of only a few, 
which would tend to tell us that using only the observed 
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Fig. 11. Left: Electron flux from observed local SNRs, with associated uncertainty band (due to observational uncertainties on 
ages, distances, radio fluxes and spectral indices). Right: Propagation effects on the electron flux originating from local SNRs, 
using either the observational constraints or the average flux and index from Fig. 2] In both panels, we have assumed a source 
cut-off energy of 10 TeV. 




Fig. 12. Primary electron flux from local sources located within 2 kpc from the Earth. Left: from know SNRs (radio constraints). 
Middle: from unknown SNRs that should accompany unassociated local pulsars (generic spectrum for each, age and distance from 
the pulsar companion). Right: from known pulsars (generic spectrum, observational constraints on age and distance). Here, all 
generic spectra refer to a spectral index of 2 and an energy cut-off of 3 TeV. 



sources translates into a slight underestimate of the actual 
local electron budget. This makes sense, since observations 
favor the brightest objects, and also since rapidly-rotating 
magnetized neutron stars with magnetic axes pointing away 
from the Earth cannot be observed as pulsars. Our samples 
of local SNRs and pulsars are shown in the age-distance 
plane in Fig. O 

In Fig. [TOl we show the electron flux obtained in the 
med propagation setup, using the local SNR properties 
summarized in Table IC.ll and assuming different cut-off en- 
ergies. The first important comment to make is that the 
whole flux is far from being described simply by a smoothly 
injected power law, since many spectral wave-like features 
are evident. Moreover, because of the interplay between the 
age and the maximal energy (see Sect. 12. 2p . we see that 
varying the energy cut-ofF from 1 to 100 TeV, though the 



latter value is probably not realistic and too high, has con- 
siderable effects; not only do new contributions arise at high 
energy when the cut-ofF value increases, but the hierarchy 
among other sources is also altered. This illustrates an addi- 
tional source of theoretical uncertainty, beside those we dis- 
cussed in Sect. 14.31 Figure [TT] illustrates the theoretical and 
observational uncertainties affecting both the source model- 
ing and the propagation modeling. The left panel allows us 
to quantify the impact of the observational uncertainties on 
the ages, distances, radio fluxes and spectral indices. This 
plot portrays the results of 1,000 Monte Carlo realizations 
in which we drew each parameter according to a flat dis- 
tribution within the observational errors. The right panel 
exhibits (i) the impact of varying the propagation parame- 
ters on the overall local SNR contribution and (ii) the differ- 
ences associated with different injection spectrum prcscrip- 
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tions. For the latter point, we compared the flux obtained 
with generic spectral properties, namely a spectral index of 
2 and a fixed normalization of Qq = 3.9 x lO^^GeV"^ (see 
end of Sect. 14. 1)) , with the observationally constrained pre- 
dictions. We again see that the global spectral shape is far 
more complex than a mere power law, and that the over- 
all flux can vary within a factor of 2-5 depending on the 
energy. 

We plot the results obtained with a template calcula- 
tion for all local electron primaries in Fig. 1121 where we 
used the med propagation setup with the Ml ISRF model 
(the large difference between the full relativistic calcula- 
tion and the Thomson approximation is evident in the 
plots). The three panels from left to right show the con- 
tributions of local known SNRs, from non-observed SNRs 
associated with observed pulsars, and from observed pul- 
sars, respectively. For the non-observed SNRs, we consid- 
ered the distances and ages of the associated pulsars, and 
we assumed an injection spectrum with an index of 2 and 
a non-observable radio flux of B{1 GHz) = 1 Jy. While the 
contribution from the non-observed SNRs is shown to be 
negligible with respect to the two others, it is interesting to 
note that our local pulsar modeling leads to a higher pri- 
mary electron flux than local SNRs. For pulsars, we sup- 
posed than 10% of the spin-down energy was converted 
into electron-positron pairs when defining the individual 
normalizations (see Sect. [S]). Although the injection mech- 
anism is subject to large theoretical uncertainties, it is still 
rather surprising to find that even when accounting for ob- 
servational constraints for pulsars, their local population 
can contribute as many primary electrons as known local 
SNRs. However, these calculations are subject to very large 
theoretical uncertainties, so no strong conclusions should 
be drawn: by no means should they be considered as pre- 
dictions, only trends. The main conclusion at this stage is 
that nearby sources dominate the flux above energies of few 
tens of GeV, the spectral imprints of which are very diffi- 
cult to predict, because of large theoretical uncertainties 
in their modeling. Finally, it is clear from the plots that 
the Thomson approximation can lead to large errors in the 
predictions. 

4.4.3. Electron flux: sum of distant plus local sources 

We can now derive a full calculation including all known 
local sources within a distance of 2 kpc from the Earth and 
a more distant smooth component. For the latter, we use 
the L04 spatial distribution, but in contrast to what we 
did in Sect. 14.4.11 we have to apply a radial cut-off to this 
smooth contribution within 2 kpc to the Earth, to avoid a 
possible double counting of the local sources. This means 
that the distant component has more impact at low energy 
than at high energy. 

Furthermore, we can also compare our results with the 
available observational data about electrons. The top left 
panel of Fig. [T3] shows an example of a calculation includ- 
ing all the components discussed above in a self-consistent 
manner, using the med propagation setup. To proceed, we 
fixed the entire parameter set from available observational 
constraints, and we tuned the other parameters to provide 
rough agreement between the calculated flux and the data. 
In particular, we adjusted the global normalizations and 
spectral indices of the two distant smooth components (one 



for SNRs, another for pulsars) and of the local pulsars. The 
parameters that we used are summarized in Table [5] 



5. Primary positrons 

The increase in the positron fraction inferred from 
the PAMELA data is barely consistent with what is 
expec ted for secondary positrons {e.g. iDelahave et al.l 
I2009D . except for th ose secondaries produ c ed and ac- 
celerated in sour c es dBerezhko et al. 



iMertsch fc Sarkad 120091: lAhlers et al 



200l iBlasil [20091: 
20091) . Additional 



spatia l arguments might still be helpful in dealing with this 
issue (jShaviv et al.ll2009f ). 

Although exotic primary contributions from dark 
matter annihilation could contribute to the positron 
budget, the most popular dark matter particle candi- 
dates can hardly exceed the secondary background un- 
less the annihilation rat e is boosted su bstantiall y (e.g. 
iBahz fc Edsio 1998: L a valle et al.l l2008at [Pelaha ve et al.l 
boost iBerestrom et al.lfeOOSi ) or the expansion rate is more 
rapid than expected in the e arly universe , before the pri- 
mordial nucle osynthesis (se e Salatil |2003| for the original 
idea, and e.g. ICatena et al.l [20091 more speciflcally for the 
positron channel) . Furthermore, the potential enhancement 
provided by dark matt er substructures was de monstrated 
to b e too small (see iLavalle et al.l l2008bl [al: iPieri et all 
[2009I) . and any other type of global enhancement was shown 
to be severely restricted by the companion upp er limit 
to the antiproton yield {e.g. iDonato et al.l l2009t ). which 
fully applies to the most popular dark matter models that 
have in no case exclusive couplings to leptons. We there- 
fore emphasize that generic dark matter candidates are 
not expected, unfortunately, to manifest themselves in the 
positron spectrum in the GeV-TcV energy range. 

For other potential sources of primary positrons, we 
note that since the seminal work on elect ron-pos i tron p air 
production in strong magnetic flclds by lErberl (jl966l ). a 
particular class of cosmic-ray sources has long been pre- 
dicte d to pr ovide electron-positron pairs: pulsars (jSturrockl 
119701 1197111 . Interestingly enough, a detailed discussion 
about the increase in the measured positron fraction and 
the potential r ole of pu l sars w as already conducted two 
decades ago bv iBoularesI (|l989lPI . Subseq uent deeper stud- 
ies have been performed since then {e.g. lAharonian et all 
119951: lAtovanet all 119951: IChi et al.| [l 996l). which were re- 
cently revisited by several authors (e.g. 'Ho oper et al] 20091 : 
iProfumol 120081: lYuksel et al. .2009; Malvshev et al.ll2009[ ). 
Since pulsars are commonly observed in the vicinity of the 
Earth, it is very likely that these astrophysical sources of 
positrons are the origin of most of the local cosmic ray 
positrons. Nevertheless, many uncertainties remain in the 
characterization of pulsars, and there has been little theo- 
retical progress in the past decade. 

Pulsars are rapidly-rotating magnetized neutron stars, 
the rotation axis of which is misaligned with the axis of 
the magnetic field, which produces a pulsed emission for 
an observer located on the cone scanned by the magnetic 
axis. They are potential end products of CCSN explosions, 
such as black holes for the most massive progenitors and 
other neutron stars unseen as pulsars (if the observer is 
always off-magnetic axis, the rotation and magnetic axes 



* iBoulare^ (| 19891 ) did also mention that dark matter was an- 
other, though exotic, possibility. 
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are aligned, or the neutron star is either non-rotating or 
non-magnetized). All pulsars must in principle be associ- 
ated with a companion SNR. This is supported by ob- 
servations since pulsar wind nebulae (PWNe) are some- 
times found inside or close to the shells of SNRs (see e.g. 
iGaensler fc Slanell2006l:lBucciantinill2008l for recent reviews 
on PWNe), though many pulsars remain non-associatec^. 
Thus, a more general picture of pulsar modeling, at least for 
cosmic-ray electrons and positrons, should involve a SNR 
as the main energy supply, and a pulsar and its PWN 
as a subsystem inject ing additional high energy electrons 
and positrons (see e.g. iBlondin et aLll2001[ ). Note that even 
though a SNR should be associated with each pulsar, the 
reverse is not true: explosions of SNela leave no compact 
objects at all, and some CCSNe do not become pulsars. 
Likewise, we stress that the cosmic-ray acceleration mecha- 
nisms are different in SNRs and in PWNe, leading to differ- 
ent observational properties. In the former case, accelera- 
tion takes place at non-relativistic shocks, while at relativis- 
tic shocks in the latter case, which might imply different 
spectral shapes for the accelerated electrons. Moreover, pul- 
sars convert a significant part of their spin-down energy into 
electron-positron pairs that are accelerated by the PWNe; 
this does not occur inside SNRs, where the energy fraction 
transferred to electrons is much lower (see Sect. 14. 

In this section, we complete our calculations of the elec- 
tron and positron flux at the Earth by including the con- 
tribution of pulsars. We adopt the same methodology as 
for SNRs, and consider two populations, (i) the local one, 
for which we constrain the individual properties from ob- 
servational data, and (ii) a more distant and smoothly de- 
scribed population. As in most recent studies, the individ- 
ual p roperties of pulsars are derived from the ATNF cat- 
alog (jManchester et aI1l2005l ). In Sect. 15.11 we sketch the 
generic model that we have adopted. In Sect. 15.21 we discuss 
the contributions of both populations to the local electron 
and positron flux. We then compare some template calcula- 
tions to the current measurements of the positron fraction 
in Sect. 



5.1. Generic pulsar modeling 

Independently of the sp ecific pulsar model, such as fo r 
instance the p olar gap (iRudernian fc Sutherland! I1975D . 
the outer gap JCheng et al.l 1976L 19861) . or the slot gap 
(jHarding et al.l 120081) models, we may summarize the 
physics relevant to the production of cosmic-ray electrons 
and positrons as follows. Gamma rays can be generated in 
the pulsar magnetosphere from inverse Compton interac- 
tions of electrons accelerated along the strong and rotating 
magnetic field with local synchrotron radiation, which can 
in turn produce electron-positron pairs by annihilating with 
photons from the local radiation fields. Those gamma rays 
can be observed as pulsed emissi on, such as those r ecently 
detected with the Fermi satellite ()Abdo et al.ll2009"al ). which 
may therefore be used to constrain the pair production. 
These electron-positron pairs are then accelerated within 
the surrounding and expanding shocked medium, namely 



^ Pulsars are likely usually expelled from the SNR system, 
but the relative velocity between both objects < 10* km/s in- 
duces a distance < 100 pc between them by the SNR lifetime. 
This distance may be larger for old systems, which makes clear 
associations more difficult. 



the PWN, located inside or offset from a more extended 
SNR. Observations of young systems such as the Crab 
nebula tell us that this acceleration can be very efficient 
and lead to huge Lorentz factors, up to '-^ 10® — 10^ 
([Atovan fc Aharonianlll996l ) . What is important when try- 
ing to predict the electron-positron yield from a pulsar is 
not their energy distribution and density close to the mag- 
netosphere, but instead their final characteristics after ac- 
celeration has proceeded and when the particles are released 
int o the ISM. This has a lready been noted and described 
bv lMalvshev et all ()2009t ). It is therefore rather difficult to 
provide accurate predictions when disregarding the whole 
dynamics at stake there, and, in this part, we mostly aim 
to survey the roles of the main ingredients that charac- 
terize pulsars rather than making peremptory predictions. 
Indeed, we show in the following that current uncertain- 
ties still make it difficult to derive anything but qualitative 

predictions. 

Fo llowing the arguments developed in iMalvshev et al.l 
(|2009[ ). to which we refer the reader for further details, we 
define the source term associated with any single pulsar to 
be 



Qp{E, X, t) = qp{E, t^) Sit - iU + SU)) Six - X,) 



(56) 



where t^, and x^, are the pulsar age and position, respec- 
tively, and Sti, accounts for a certain delay in the release 
of cosmic ray electrons into the ISM after the supernova 
explosion. Generic pulsar s should have ceased th eir PWN 
phases after ~10-100 kyr (jGaensler fc Slan3l2006|) : for sim- 
plicity, however, since characterizing the PWN evolution 
is far beyond our purpose here, we assume that St^, = 
in the following. T he pulsar age is usually estimated from 
the spin- down age (jOstriker fc Guiinlll969[ ). which is only 
relevant to the spin-down magnetic radiation approxima- 
tion, that involves the rotation period P and its first time 
derivative P, given by 



t 



pulsar 



p 

2P 



(57) 



As in the SNR case, we again emphasize that this age esti- 
mate relies on current observations, so that the actual age 
used for cosmic-ray transport calculations should have an 
additional d/c term, where d is the distance of the pulsar to 
the observer. We stress, however, that using the spin-down 
age for the pulsar age estimate turns out, in many cases, to 
be erroneous (see e.^.^Gacnsler & Frail 2000). Nevertheless, 
for simplicity, we adopt this method to deal with local pul- 
sars in the following. 

For the energy spectrum, we adopt the same general 
shape as used previously for SNRs (sec Ea. H^ . i.e. a power- 
law of index 7 with an exponential cut-off at energy E,, as 



^" 'I 



I -e: 



(58) 



In contrast to the SNR case for which the spectral index 
can be constrained from radio observations, the spectral 
index associated with high energy electrons from pulsar 
can hardly be constrained by the radio observations of the 
pulsed emission. Indeed, this pulsed emission originates in 
regions close to the pulsar magnetosphere, where the ac- 
celeration processes are not yet achieved. An alternative is 
to use the spectral indices derived from PWN observations. 
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when available. To simplify the discussion, we use 7 = 2 in 
the following, unless other values are specified. 

The normalization Qo is intimately linked to the total 
rotational energy Wq of the pulsar, a fraction / of which is 
released in the form of electron-positron pairs, such that 

/•OO 

/ dEEqp{E) = fWo, (59) 

where Wq can be constrained from measurements assuming 
that the whole energy lost is carried by the magnetic dipole 
radiation, such that 

Wo^ETd,Jl + —) , (60) 

V TdccJ 

where E is the spin-down luminosity and Tdoc = Eq/Eq 
is the typical pulsar decay time. Note that the index v 
indicating the age dependence is in principle related to the 
braking index k that defines the rotation deceleration f] 
— fi'^. where fi is the angular velocity, in terms v ~ (k + 
l)/(fc — 1). In the spin-down approximation, k = 3, and 
therefore v = 2. Nevertheless, k can also be computed if 
the second time derivative of fl is known, k = —flQ/fl^. In 
that case, it is usually found to slightly differ from 3. 

This illustrates again the large degree of theoretical un- 
certainties arising when trying to model pulsars, even in 
simple approaches. We still use the spin-down approxima- 
tion, and therefore assume that v = 2 in the following. 

Since accounting for the details in the pulsar model- 
ing is beyond the scope of this paper, we adopt the source 
model defined by Eqs. ((551 - [50]) . and assume a universal 
decay time of Tdoc = 1 kyr. Using this latter input and 
the ATNF data for th e rotation period and its derivative, 
iMalvshev et~all (|2009D found typical values of Wq - 10"^ 
erg, i.e. one or two orders of magnitude below the char- 
acteristic supernova energy release, in agreement with the 
picture of a pulsar as a subdominant energy supply beside 
its companion SNR. Finally, it is usually supposed that a 
rather significant fraction of the spin-down energy is con- 
verted into electron-positron pairs that are accelerated in 
the PWNe. In the following, we assume by default that 
this fraction is / = 0.1, unless otherwise specified. Though 
the typical spin-down energy of pulsars is a hundred times 
lower than the typical kinetic energy of SNRs, the rather 
efficient pair conversion occurring in pulsars implies that 
the overall energy carried by charged leptons is similar in 
both cases — we recall that the fraction of kinetic energy 
converted into leptons is around 10^^ — 10^^ in SNRs (see 
Sect.O]). 

For the nearby known pulsars, we extract the distances 
and the ages from the ATNF catalog. For more distant 
objects, say above a few kpc, we can safely use the contin- 
uous limit as we did for SNRs, and write the corresponding 
source term as 

QpiE,x)^qpiE)f{x) . (61) 

The energy spectrum is assumed to be the same as in 
Eq. f|58p . and the normahzation Qo is defined by an ex- 
plosion rate that is similar to the ease of SNRs. 

5.2. Local versus distant pulsars 
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Fig. 13. Selection of local ATNF pulsars (within 2 kpc) in the 
P ~ P plane. Only those pulsars with a spin-down age lower 
than 30 Myr were selected in our sample. 



Wc selected local pulsars from the ATNF catalog, im- 
posing a few constraints. First, we applied a radial cut-off 
of 2 kpc, as for the SNR treatment, farther pulsars being 
accounted for with a smooth spatial distribution modeling. 
For the latter, we adopted the L04 model, as for distant 
SNRs (see Sect. I4.4TT] for more details). Considering ener- 
gies above the GeV scale imposes an upper limit of the 
pulsar age < 30 Myr, which decreases to 1 Myr above 100 
GeV (see Fig. [7]). We therefore restrict our sample by re- 
quiring that the ages < 30 Myr. This selection procedure 
is depicted in Fig. [131 where the local ATNF pulsars are 
reported in the P-P plane. Our final sample contains a bit 
more than 200 objects located at fewer than 2 kpc from 
the Earth, which is quite large compared to the number of 
observed SNRs, i.e. 27, discussed in Sect. 14. 4T^ 

The positron flux derived from this pulsar selection is 
obviously identical to the pulsar contribution to the elec- 
tron flux discussed in Sect. 14.4.21 because of the pair pro- 
duction mechanism, and is shown in the right panel of 
Fig. [T2j We remind the reader that a SNR counterpart was 
systematically added to each local pulsar, except for pulsars 
that already had an observed and identified SNR counter- 
part featuring in our SNR list. In case of a non-observed 
SNR counterpart, we assumed a radio fiux B{\ GHz) = IJy 
to set its spectrum normalization (see Sect.lllj). These non- 
observed SNR counterparts contribute only to the electron 
fiux, which is reported in the middle panel of Fig. 1121 

As briefiy mentioned in Sect. 14.4^ using the spin-down 
approximation to constrain the energy released in the form 
of electron-positron pairs leads pulsars to be important 
sources of local high energy electrons and positrons, as in- 
tense as the observed local SNRs. We employed a pair con- 
version efficiency of / = 0.1 in this calculation, which might 
be optimistic, but still, decreasing this efficiency to a few 
percent would lead to a quite significant contribution to the 
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Table 4. Main positron sources among the ATNF nearby pulsars. We rank the pulsars from the largest contribution to the flux 
in different energy bins, assuming a spectral index of 7 = 2 and a cut-off energy of 1 TeV. All other parameters are derived from 
the ATNF catalog. 



local electron and positron budget. It is noteworthy that the 
local positron flux is dominated by a few objects among our 
^--^200 selected objects. The main sources and their prop- 
erties are listed in Table |31 Although pulsar modeling is 
subject to many and large theoretical uncertainties, and 
despite the simplistic model we have employed to set the 
individual normalizations, our results suggest that SNRs 
might not be the only prominent sources of electron and 
positron cosmic rays. Nevertheless, this complete approach 
involving pulsars and SNRs in a self-consistent framework 
should be studied in greater detail, which may be promis- 
ing for gaining deeper insights into the understanding of 
cosmic ray leptons. 

A full template calculation of the positron flux in- 
cluding secondaries and primaries from distant and lo- 
cal pulsars is compared with existing data in the bot- 
tom left panel of Fig. [HI We used t he data from HEAT 
(|DuVernois et all l200llPcAPRICE (iBoezio et al.l 1200011 . 
and AMS-01 (jAlcaraz et al.l |2000D . which provide con- 
straints at energies lower than ~20 GeV only. From this 
plot, we see that secondaries can already account for a large 
fraction of the low energy positron flux, which constrains 
the properties of the additional primaries. We have in fact 
some freedom to tune the normalization and the spectral 
indices of the distant and local component, since they are 
poorly constrained. For the normalization, we have to ad- 
just the pair conversion efficiency / to values lower than 
0.1, which reinforces the idea that pulsars can contribute 
to the local electron and positron flux quite naturally. This 
is discussed further when tackling the positron fraction case 
in Sect. Likewise, there is no particular need to invoke 
hard spectral indices < 2. The parameters that we have 
used are summarized in Table [5l 

Finally, as for local contributions to the electron flux, we 
again point out that using the Thomson approximation to 
deal with the energy losses can lead to large errors and fake 
predictions of very peaky features in the overall positron 
spectrum. 



6. Full electron and positron results — discussion 

In this section, we now perform a template calculation with 
reasonable parameters for all previously discussed ingredi- 
ents to demonstrate that pure astrophysical processes easily 
account for current measurements. It is quite simple to find 



configurations in good agreement with the data of the sum 
of electrons plus positrons or/and on the positron fraction, 
but it is far more difficult to find additional agreement with 
the separate electron and positron data. Many other studies 
have focused only on a few of these data sets, more rarely 
on all together. 

As an important and preliminary remark, we emphasize 
that playing with the parameters associated with very few 
local sources within reasonable ranges allows us to fit all 
the data quite easily, though roughly, simply by adjusting 
the low/high energy source hierarchy. Nevertheless, for the 
sake of illustration, we attempt to use the full set of ob- 
servational constraints that we dispose of without turning 
around. This is to illustrate how difficult and complex this 
self-consistent exercise can be. However, to simplify the dis- 
cussion, we stick to the mtd set of propagation parameters, 
keeping in mind that uncertainties affecting propagation 
lead to 1 order of magnitude uncertainty in the overall 
predictions. For the solar modulation, we use a Fisk poten- 
tial of 600 MV. 

In addition to our template calculation, the data them- 
selves, when assumed trustworthy, can provide the actual 
separate positron and electron fluxes at the Earth. If we 
define /_ and /+ as the electron flux and the positron flux, 
respectively, and assume that some fitting functions /tot 
and /frac for the measured total flux {e.g. from Fermi data) 
and for the positron fraction {e.g. from PAMELA data), 
respectively, exist and are known, we then readily find 



f+{E) = fi^^E) X f,ot{E) 
f^{E) = hoi{E){l- ff^UE)) 



(62) 
(63) 



Although we do not perform the exercise in the following, 
we point out that obtaining individual predictions in good 
agreement with /+ and /_ in the relevant energy range 
automatically ensures close agreement with the measured 
positron fraction and total electron plus positron flux. This 
might sound tautological, but is still a helpful method if 
one wishes to optimize her/his preferred CR model. 

6.1. Electron and positron spectra: separate fluxes and sum 

We compare a template calculation of the electron and 
positron fluxes with the available data. For electrons and 
positrons, as measured separately, we use the data from 
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Table 5. Injected energy, converted fraction, spectral indices, and cut-off energies used for the overall template electron and 
positron flux calculation. The symbol f indicates that we used observational constraints. For local SNRs, we used a global extra- 
factor of 6, which corresponds to assuming a magnetic field of ~ 30 instead of 100 /xG in Eq. (|50p . 



CAP RICE (iBoezio et al.l [200?ih. HEAT (|DuVernois et al.l 
I2001D . and AMS-01 (jAlcaraz et al.ll2000( ). Note that these 
data are quite poorly constraining at high energy, since 
they refer to energies below a very few tens of GcV, and are 
affected by large statistical errors for the highest available 
energy bins. Future data from PAMELA and AMS-02 may 
certainly provide much tighter separate constraints on 
both components. For the sum of elec trons and positrons , 
we employ the data from Fermi (lAbdo et al.1 l2009bt 
iPesce-Rollins fc for the Fermi-LAT Collaboration! I200E " 



ATIC (IChang et al.l I2008D and HESS ()Aharonian et alj 
120081 12009air ir is meaningless to perform a fit, given the 
large theoretical uncertainties and the large number of 
parameters. We therefore empirically select the parameter 
sets that provide the smallest statistical errors whenever 
possible. 

The top and bottom left panels of Fig. [14] show tem- 
plate results for the electron flux and the positron flux at 
the Earth, respectively. The parameters that we used to 
calculate all contributions are reported in Table \E[ 

For electrons, we see that the data of CAPRICE and 
AMS-01 are barely consistent, but we can identify a me- 
dian configuration. Our template model is partly typified 
by a smooth SNR-plus-pulsar component (following the L04 
spatial distribution) that fits the low energy part — the 
smooth pulsar contribution is actually completely negligi- 
ble. We used a rather soft spectral index of 2.4 to shape 
this low energy smooth contribution, and a rather large 
value of r^; (still consistent with theoretical or observational 
constraints — see Sect. 14. ip . Nevertheless, we point out 
that using another spatial distribution, for instance with a 
stronger /softer spatial gradient, may have given the same 
result with other appropriate spectral indices and normal- 
izations (see right panel of Fig. [5]). This smooth component 
rapidly declines for energies higher than few tens GcV be- 
cause of the 2 kpc radial cut-off around the Earth that 
we imposed before adding the local component. Above 10 
GeV, we see that the local SNR contribution comes into 
play, clearly dominating over local pulsars. We underline 
that this hierarchy can hardly be considered as a generic 
prediction, since both types of sources are expected to pro- 
vide similar amounts of electrons (see left and right panels 
of Fig. [12]). Moreover, we adjusted the local SNR flux by a 
global factor of 6 with respect to our constraint recipe given 
in Eq. (|50|) . which corresponds to using a mean magnetic 
field of ^ 30 fiG instead of 100. We remark that only a few 
sources dominate the overall flux (see Fig. [T^ and Table[3]), 
so modifying the parameters of these latter only may have 
a very strong impact on the predictions. We also note that 
another pulsar to SNR contribution ratio could easily lead 
to a similar overall shape for the electron flux. Likewise, 



slightly relaxing our chosen constraints either directly on 
the individual intensities or on the ages and distances, even 
across quite limited ranges, would allow far more freedom to 
shape each contribution differently. Nevertheless, the source 
hierarchy is here imposed by the observational constraints, 
which are relevant to Eq. ([50]) in the SNR case and Eq. (jBU]) 
in the pulsar case; other methods based on different as- 
sumptions and different data would obviously provide dif- 
ferent results. Finally, we stress that the local source flux is 
not characterized by individual prominent spiky peaks, as 
would be the case in the Thomson approximation for the en- 
ergy losses (see Fig. IT^ . Additional constraints on sources 
of electrons come from the sum of electrons and positrons 
and from the positron fraction, as we later discuss. 

For the positron flux, the important point is that sec- 
ondaries are most likely predominant below 5 GeV. The 
corresponding calculation is less affected by theoretical 
uncertainties, so we could assume it to be a prediction. 
Indeed, the progenitor proton flux is not expected to vary 
significantly within a few kpc, and the main uncertainty 
apart from propagat ion is the average ISM gas density 
(jPelahave et al.l l2G09D . such that the low energy positron 
flux may vary by at most a few tens of %. Beyond 5-10 GcV, 
there are poor constraints on the positron flux, so that the 
contribution of primaries may have any shape. The positron 
fraction provides additional information and constraints, as 
we later demonstrate. Nevertheless, we see that our choice 
of parameters for pulsars (distant and local) makes their 
contribution significant above ~ 5 GeV. From Table [S] it 
appears that assuming less than 1% of magnetic energy 
converted into electron-positron pairs is enough to make 
the pulsar contribution arise around 10 GeV; this is rather 
independent of the injected spectrum (we here assumed a 
spectral index of 2). We also remark that the smooth pul- 
sar contribution is found to be negligible with respect to 
the local one. Finally, we emphasize that the overall spec- 
tral shape obtained for the positron flux weakly depends, 
in fact, on the assumed spectral index for pulsars. Indeed, 
as already shown in the right panel of Fig. [T21 decreas- 
ing the intensity of a few old sources (such as Gcminga or 
B1929-f-10) could easily result in a much harder spectrum 
at high energy and even a lower contribution below 10 GeV. 
Therefore, more tightly constraining the amplitude of the 
observed sources is a priority for future works. 

The top right panel of Fig. IHlshows the sum of electrons 
and positrons, consistently with the separate results that 
we have just discussed above. The most distinctive infor- 
mation one can extract from this plot is probably the energy 
cut-off that explicitly appears around a few TeV, which ob- 
viously encourages us to set a cut-off energy around in this 
range. In our template calculation, the high energy part 



Delahaye et al.: Galactic electrons and positrons at the Earth 



25 





Fig. 14. Template calculation including all primary (discrete local and smooth distant) and secondary electrons and positrons in 
a self-consistent modeling, using the med propagation setup. Top left: electron flux. Top right: electron and positron flux. Bottom 
left: positron flux. Bottom right: positron fraction. 



of the spectrum, above ^ 50 GeV, is actually dominated 
by our local SNR sample. This cannot be considered as a 
robust prediction, since we could have obtained a dominat- 
ing pulsar contribution by slightly modifying the injection 
parameters. For instance, we could have reduced the local 
SNR yield, both depleted the supply from Geminga and in- 
creased that from Vela, and finally allowed a larger fraction 
of ^ 2% of pair conversion. 

We point out that the current constraints on sources are 
far too weak to ascertain the predictive power of our tem- 
plate calculation. Nevertheless, it is clear that the source 
modeling is the key point in the understanding of the high 
energy CR electrons, at least for identifying more clearly 



those few sources which may set most of the local flux. 
We also underline that it is important to challenge any CR 
electron prediction with the separate electron and positron 
data, since reproducing the sum does not necessarily en- 
sure the relevance of the model. It will be important in the 
future to have independent sets of far higher quality data; 
PAMELA and AMS-02 carry many hopes. 

6.2. Positron fraction 

Measurements of the positron fraction contribute ad- 
ditional constraints that should be fulfilled consis- 
tently with the previous data. In the bottom right 
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panel of Fi^. \T^ we h a ve used th e dat a from 
HEAT tearwic k etall Il997t iBeattv et all |2004[ ). AMS- 
OljAMS-Ol C ollaboration et all l2007ir ~and PAMELA 



(jPAMELA Collaboration et al.ll2010D . From this plot, our 
template calculation turns out to be consistent with the 
data, especially above 10 GeV. At lower energies, there 
are large discrepancies among the available measurements, 
which renders the interpretation complicated, though a 
more realistic treatment of the solar modulation — which is 
beyond the scope of this paper — might improve their level 
of consistency. Note that we can accommodate a slightly 
increasing positron fraction with very reasonable parame- 
ters for the local pulsar modeling, i.e. quite modest val- 
ues of the converted energy fraction in e"*"e~ pairs and 
spectral indices, without the need to tune the parameters 
for the sources individually. Indeed, for our local samples 
of SNRs and pulsars, we emphasize again that we con- 
strained the whole contributions from observational con- 
straints and fixed assumptions. It would have been quite 
easy to make this fraction increasing much rapidly with en- 
ergy by considering additional assumptions for the pulsar 
power. We instead considered the simplest approach of the 
spin-down approximation for which the injection rate scales 
like the squared age (see Eq. [nO])- Had we adopted alterna- 
tive modeling with a shallower dependence on the age (sec 
e.g. IChi et al . 1996: Zhang & Cheng 2001; Profumo 2008), 
we would have placed more weight on younger objects, and 
thereby to the high energy part of the spectrum. As for the 
other measurements discussed above, we conclude that the 
positron fraction can be explained by pure astrophysical 
processes, even if the predictive power of our calculation 
remains weak. 



6.2.1. Final comments 

We wish here to highlight a few important points. First, 
we can fairly state that standard astrophysics is capable 
of reproducing the existing data, despite the large theo- 
retical uncertainties associated with the fiux calculation. 
Second, we attempted to convince the reader that an over- 
all check of the separate electron and positron flux calcula- 
tions against the associated separate data is mandatory to 
determine the relevance of a model: this step required the 
use of the four sets of data used in Fig. [14] separately. This 
was not covered in the literature so far, where most authors 
focused on either the sum of CR electrons and positrons or 
on the positron fraction. Third, we have shown that the full 
relativistic treatment of the energy losses could lead to a 
global spectral shape devoid of peaky structures. Fourth, we 
have discussed a template example where we fixed the local 
source hierarchy based on a few assumptions and observa- 
tional constraints. It is clear, however, that a more refined 
analysis with a far more accurate source modeling will help 
us to reduce the theoretical uncertainties and clarify the lo- 
cal source hierarchy; it will also probably lead to different 
results. Nevertheless, we have proposed an exhaustive set 
of analytical tools together with a robust method to tackle 
such a detailed analysis. Finally, we stress that we have 
derived our template calculation using our med set of prop- 
agation parameters to facilitate the reading of the plots, 
and have therefore not discussed the uncertainties originat- 
ing in the propagation side, which are actually large, as 
discussed in Sect. 14.3] for individual sources — see bottom 
right panel of Fig. [T] and Fig. [TT] 



7. Summary and conclusions 

We attempted to perform an exhaustive study of the main 
ingredients relevant to the calculation of the high energy 
CR electron and positron fluxes in the GeV-TeV energy 
range. We have underlined the complexity of the source de- 
scription, by emphasizing the potentially strong effects of 
nearby sources, and shown that though we can hardly pre- 
dict the local electron and positron fluxes with accuracy, 
reasonable parameterizations of the source and the propa- 
gation modeling can fairly well and simultaneously account 
for the current measurements of (i) electrons, (ii) positrons, 
(iii) the sum of them, and (iv) the positron fraction. We 
summarize hereafter the way in which we have proceeded. 

In Sect. [21 we presented our propagation modeling in de- 
tail. We reviewed the analytical solutions to the transport 
equation in Sects. \TJ\ and and derived explicit links be- 
tween the source spectrum and the propagated spectrum 
in Sect. 12.31 More importantly, we developed in Sect. 12.41 
a method to account for the relativistic energy losses, pro- 
viding useful fit formulae that can be used for any inverse 
Compton processes involving a black-body photon target. 
We have indeed shown that the Thomson approximation 
was by far insufficient and could lead to fake predictions 
of peaky signatures in the electron spectrum as potential 
imprints of local sources. In Sect. 12.51 we recalled our prop- 
agation parameters with emphasis on the energy losses. In 
particular, we developed a quite tightly constrained local 
ISRF modeling and explained in detail how its different 
components result in decreasing steps in the overall energy 
loss function, as a consequence of Klein-Nishina effects. 

In Sect.[3J we revisited the predictions for the secondary 
electron and positron fluxes, finding slight differences from 
the p revious analysis we performed in [Delahave et al.] 
(]2009l ) — for positrons only, and using the Thomson ap- 
proximation for the energy losses. These differences are 
caused by our novel relativistic treatment of the energy 
losses. We found a slightly harder secondary spectrum, with 
an overall flux higher by - 10% at 100 GeV up to ~ 50% 
at 1 TeV. We emphasized, however, that this can still not 
explain an increasing positron fraction. We provided fit- 
ting formulae of our results in Appendix JB] valid in the 
med propagation setup. 

In Sect. HJ we described in detail the calculation of the 
primary electron flux. We explored thoroughly the possi- 
ble ways of normalizing the injection rate of CR electrons 
in Sect. 14.11 focusing on a smooth distribution of sources 
as well as on discrete sources, and showed that the theo- 
retical uncertainties are impressively large, reaching 2 or 3 
orders of magnitude. In particular, the fraction of energy 
transmitted to CR electrons is quite an issue. Nevertheless, 
we sketched a method to normalize an event-like source in- 
jection density from its measured radio flux, which allows 
us to use observational constraints whenever available. We 
demonstrated in Sect. 14.21 that any smooth spatial distri- 
bution of sources could be approximated with a mere z- 
exponential disk up to very high accuracy for electrons of 
energy above a few tens of GeV, showing that the spe- 
cific spatial distribution had negligible impact on the pre- 
dictions, for which only the local density of sources really 
accounts. Since local effects are more important above a 
few tens of GeV, we exhaustively studied in Sect. 14.31 the 
impacts of the parameters relevant to the local event-like 
source description, i.e. the age, the distance and the spec- 



Delahaye et al.: Galactic electrons and positrons at the Earth 



27 



tral index, as well as of the propagation parameters (see 
Fig. [71) . We notably illustrated how irrelevant the individ- 
ual spectral indices were to inferring the overall spectral 
shape produced by a population of non-identical objects. 
We have shown that the overall spectrum was mostly set 
by a hardly predictable hierarchy in the local sources. The 
energy range characterizing a single contribution is mostly 
bounded from above by the age and cut-off energy, and from 
below by the distance. For the former bound, we demon- 
strated that the Thomson approximation led to erroneous 
predictions, with significant underestimates of the maximal 
energy and fake peaky features. We also demonstrated that 
changing the diffusion coefficient had a strong impact on 
the results. Finally, in Sect. 14. 4[ we discussed a template 
calculation including the constrained contributions of all 
known local sources lying within a distance of 2 kpc from 
the Earth, showing that pulsars might contribute as many 
primary electrons as SNRs. We again emphasized the quite 
poor predictive power that can be achieved even when ac- 
counting for observational constraints. This is mostly due 
to the theoretical uncertainties affecting the source model- 
ing. Nevertheless, we indicated a few objects that probably 
dominate the overall high energy flux. 

In Sect. [51 we studied the contribution of pulsars to 
the primary positron flux, using the spin-down approxi- 
mation to constrain the injection rate of local pulsars. As 
for primary electrons, we demonstrated that the very large 
theoretical uncertainties makes it difficult to claim robust 
predictions. 

For the local sources, we proposed a novel approach in- 
cluding not only known SNRs or pulsars, but also adding 
a SNR counterpart to each non-associated pulsar. Indeed, 
since a pulsar is the relic neutron star of a core-collapse su- 
pernova explosion, it should be accompanied by a SNR for 
consistency reasons. We accounted for these non-observed 
SNRs needing to have a radio flux below the current exper- 
imental sensitivities when normalizing their injection rates 
of electrons. In contrast, all SNRs are not expected to have 
a pulsar companion, in particular those SNRs coming from 
SNla explosions. A more refined and realistic modeling of 
sources including such composite contributions and further 
accounting for the dynamics of CR injection remains to be 
considered in terms of electron flux. This is, however, far 
beyond the scope of this paper, though we attempted to 
pave the road for such an ambitious study. 

Finally, in Sect. [51 we presented a template calculation 
including all sources of primaries, in which we used obser- 
vational constraints whenever available (see Fig. [H]) . We 
considered the true complexity of the source description 
within a self-consistent framework, at variance with many 
other studies that merely added local contributions to other 
so-called standard predictions. We demonstrated that for 
reasonable assumptions about the parameters of both the 
local and distant components, all of the discussed astro- 
physical processes can account rather well for all measure- 
ments of cosmic ray electrons and positrons independently. 
We pointed out that comparing any prediction with only 
the positron fraction or the sum of electrons and positrons is 
incomplete, and does not ensure its consistency. We empha- 
size that our template calculation cannot stand for a robust 
prediction because of the very large theoretical uncertain- 
ties underlined along this paper. Our results should instead 
be considered as a proof that current measurements are 
compatible with pure astrophysical processes, so that one 



can hardly claim any prominent anomaly at the moment. 
More importantly, we hope to have shown that a standard 
model of CR electrons, standard in the sense that depar- 
tures from observational data would sign something unex- 
pected or exotic, is far from being reached, since we even do 
not control the modeling of the local environment (mostly 
the sources), which is of paramount importance here. This 
may have drastic consequences; for instance, it indicates 
how well one can control or predict the CR electron con- 
tribution to the high energy part of the Galactic diffuse 
gamma-ray emission, or to the diffuse radio emission. In 
turn, however, this diffuse emission can provide additional 
information and constraints on the astrophysical processes 
at stake. More accurate predictions may only emerge for 
observables that involve features averaged over large spa- 
tial and/or timescalcs, to guarantee protection against po- 
tentially large and hardly controllable fluctuations in the 
calculation results. This is for instance the case for high 
energy stable nuclei. 

Such conclusions may sound disappointing but they im- 
ply many interesting by-products and open anyway a very 
broad landscape for future theoretical improvements. For 
instance, we have shown that a few local sources may domi- 
nate the electron flux. If so, far more refined models of these 
sources may allow us to use electrons as independent trac- 
ers of the local diffusion coefficient. Positive measurements 
of anisotropics in the CR electron flux might also provide 
insights in the power hierarchy among local sources, pro- 
vided this hierarchy is strong; indeed, each flux energy-bin 
is most likely filled by contributions from several sources 
located at different positions in the sky. In any case, the 
connection between the local CR electrons and the physics 
of sources is probably one of the most important issue to in- 
vestigate further in future studies; low energy protons might 
also provide additional constraints on the local sources. 

Finally, for the searches for exotic signatures, e.g. from 
dark matter annihilation or decay, the large theoretical un- 
certainties associated with the astrophysical processes pre- 
vent us from deriving much stronger constraints at the mo- 
ment, unfortunately. Nevertheless, we underline that invok- 
ing any exotic contribution is absolutely unnecessary to 
understand the current data, and is thereby quite easily 
arguable. It seems therefore likely that models which are 
over-tuned to fit the positron data, aside from independent 
motivations coming from particle physics, arc of very weak 
relevance. We recall that almost all well-studied dark mat- 
ter particle candidates, e.g. in the framework of supersym- 
metric or extra-dimensional theories, are expected to man- 
ifest themselves in neither the local positron nor aritipro- 
ton spectra (e.g. Lavalle et al.ll2008at iDelahave et al.|[2008l : 



iDonato et al.ll2004D . Were dark matter annihilation really 
to enhance the positron budget, far more work would be 
necessary to prove it, by means of multi- messenger and mul- 
tiwavelength approaches. This would require great efforts to 
face the complexity of the astrophysical backgrounds and, 
more importantly, this will have to rely on self-consistent 
calculations of both the signals and backgrounds. 
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Appendix A: Analytical solutions of the spatial 
integral 

We briefly present the solutions of the spatial integral of 
the Green function convolved with a disk-like spatial dis- 
tribution for sources. We distinguish two cases: (i) a ho- 
mogeneous flat disk, which is relevant when dealing with 
secondary CR electrons and positrons (see Eq.HT]), and (ii) 
a z-exponential disk (see Eq. I5ip . relevant when dealing 
with primaries as long as the distribution of sources does 
not exhibit too strong gradients over a distance fixed by 
the half-thickness of the diffusion zone ^ L. 
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In both cases, disregarding the radial boundary condi- 
tions, the radial solution is simply given by: 



exp 



exp 

>oo 1 

'^0 27r 



27r 



1 — exp < — 



(A.l) 



max mm 



A2 



If one wishes to take into account the radial boundary, 
it is enough to add a second propagator that will cancel out 
the first one at r = i?. The 2D-propagator becomes: 



exp^-a(r.) ^^'^ MA.2) 



where the image satisfies r^™ = rs/a(rs) with the scale pa- 

2 

rameter a{r) = This is equivalent with replacing Rq by 
and A by A-^ (except for the first A^ which is in fac- 
tor). In the case of secondary production, when the source 
term is homogeneous in the disk, the integration over the 
total disk is analytical and gives: 



llT= I I rgr{r,E^r,,E,)drde 

JQ JO 



e 
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a(i?0 
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^ 1 1 A2 ) m! I ^ 
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A.l. Vertical solution for the homogeneous disk 
approximation 

We quote the result obtained by iDelahave et alj ()2009l ) for 
secondary positrons. We have: 



horn 



(A.3) 
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The latter case corresponds to the Helmholtz solution, 
while the former is the image solution, for which we have 
z^"" = 2nL + (-l)"z„,ax and z™'" = 2nL - (-l)"z„,ax. 
Throughout the paper, we have used a disk of half-thickness 
Zmax = h = 0.1 kpc. 



A. 2. Vertical solution for the z-exponential disk 
approximation 

If we consider sources that exhibit exponential vertical pro- 
files, the integral has to be performed over the complete 
diffusion zone: 



rfz,^,(A,z = 0^z,)e(-l-l/-") 



(A.4) 
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The first case corresponds to the image solution and 6„ 
represents + (— The second case corresponds to 
the Helmholtz solution. 

Appendix B: Fitting formulae for secondary 
electrons and positrons 

Here, we provide useful fitting formulae for our template 
computations of the interstellar secondary fiux, for both 
CR electrons and positrons, in the med propagation setup 
and with the Ml ISRF model for the full relativistic energy 
losses (see sections 12.51 and [51 for more details). These for- 
mulae do not account for the solar modulation that one has 
to implement oneself, and are relevant from ^ 1 GeV up to 
a few TeV. We used the following parameterization for the 
flux, expressed in units of (GeVcm^ ssr)"-'^: 



N 



0soc(-E') = exp < 



. 4=0 



In 



E 



IGeV 



(B.l) 



For secondary positrons, we find 

af = I - 5.46298,-3.52896,-0.0887432,0.0319396, 

-0.00339393, 0.000107393| , (B.2) 
while for electrons, we have 

af = I - 6.00407,-3.40715,-0.0977409,0.0346854, 

-0.00441225, 0.00019| . (B.3) 

Note that rescaling these formulae by a global factor mostly 
means rescaling the averaged local value of the gas density 
in the disk. 

Appendix C: Local supernova remnants 

SNRs are very numerous but not always easy obj ects to 
detect . The most complete availabl e catalog is th at of lGreeij 
(|2005[ ). which was later updated (iGreenI 120091) . Moreover 
the Canadian Galactic Plane Survey "(iTavlor et all l2003l 
CGPS hereafter) has focused on many of these objects. In 
the following, we discuss the SNRs that lie less than ~ 2 kpc 
away from the Earth. We classified them by distance to the 
Earth. All the important data we have used are summed up 
in Table IC.ll Considering the number of recent discoveries 
in this region, it is more than probable that other remnants 
will be discovered in the future. 
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# 


SNR 
G+long+lat 


other name 


distance 
[kpc] 


radio index 


Brightness 

[Jy] 


age 
[kyr] 


Pulsar 


1 


18.95-1.1 




2. ± 0.1 


0.28 


40 


11.75 ± 0.85 


? 


2 


65.3+5.7 




0.9 + 0.1 


0.58 + 0.07 


52 


26+1 




3 


65.7+1.2 


DA 495 


1.0 + 0.4 


0.45 + 0.1 


5 


16.75 + 3.25 


unknown 


4 


69.0+2.7 


CTB 80 


2.0 + 0.1 


0.20 + 0.10 


60 ± 10 


20+1 


J1952+3252 


5 


74.0-8.5 


Cygnus Loop 


o.54+;i-,v« 

— 0.08 


0.4 + 0.06 


175 ± 30 


10+1 




6 


78.2+2.1 




1.5 + 0.1 


0.75 + 0.03 


275 ± 25 


7+1 




7 


82.2+5.3 


W63 


2.3 + 1.0 


0.36 + 0.08 


105 ± 10 


20.1 + 6.6 




8 


89.0+4.7 


HB 21 


1.7 ± 0.5 


0.27 + 0.07 


200 ± 15 


5.60 ± 0.28 




9 


93.7-0.2 


CTB 104A or DA 551 


1.5 + 0.2 


0.52 + 0.12 


42+7 


50 ± 20 




10 


114.3+0.3 




0.7 


0.49 + 0.25 


6.4 + 1.4 


7.7 ± 0.1 




11 


116.5+1.1 




1.6 


0.16 + 0.11 


10.9 + 1.2 


20 + 5 


B2334+61 ? 


12 


116.9+0.2 


CTB 1 


1.6 


0.33 + 0.13 


6.4 + 1.4 


20 + 5 


B2334+61 ? 


13 


119.5+10.2 


CTA 1 


1.4 + 0.3 


0.57 + 0.06 


42.5 + 2.5 


10 + 5 


JOOlO+7309 


14 


127.1+0.5 


R5 


l.+O.l 


0.43 + 0.1 


12 + 1 


25 + 5 




15 


156 2+5 7 




0.8 + 0.5 


2 0+i-i 

^•'-'-0.7 


4.2 + 0.1 


10+1 


B0450+55 


16 


160.9+2.6 


HB 9 


0.8 + 0.4 


0.48 + 0.03 


~ 75 


5.5 ± 1.5 


B0458+46 


17 


180.0-1.7 


S147 


1.2 ± 0.4 


0.75 


74 ± 12 


600 ± 10 


J0538+2817 


18 


184.6-5.8 


Crab nebula 
or 3C144 or SN1054 


2.0 ± 0.5 


0.3 


1,040 


7.5 * 


B0521+31 


19 


189.1+3.0 


IC 443 


1.5 + 0.1 


0.36 + 0.04 


160 ± 5 


30 or 4 




20 


203.0+12.0 


Nlonogem ring 


U.^oo_Q Q27 






86 ± 1 


B0656+14 


21 


205.5+0.5 


Monoceros Nebula 


1.63 + 0.25 


0.66 ± 0.2 


156.1 ± 19.9 


29+1 




22 


263 9-3 3 


Vela(XYZ) 


295 ± 075 


Vciri3,blG 


2 000 ± 700 


11.2 + 0.1 


B0833-45 


23 


266.2-1.2 


RX J0852. 0-4622 or Vela Jr 
or SN1300 


0.75 + 0.01 






3.5 + 0.8 *? 


J0855-4644 ? 


24 


276.5+19.0 


Antlia 


0.2 + 0.14 






> 1,000 


B0950+08 


25 


315.1+2.7 




1.7 + 0.8 


0.7 




50 + 10 


J1423-56 


26 


330.0+15.0 


Lupus Loop 


1.2 + 0.3 






50 ± 10 


B1507-44 ? 


27 


347.3-0.5 


SN393 


1.+ 0.3 






4.9 * 





Table C.l. Characteristics of nearby SNRs. Spectral index and brightness are inferred from measurements made at 1 GHz. 
Uncertainties in bold are not taken from bibliographic references, but just correspond to a rough uncertainty in the last relevant 
digit; hence they can be underestimated. An age is flagged with a ★ for an historical remnant; in this case, the age uncertainty is 
set from the distance uncertainty. Note that these ages are the observed ages, which differ from the actual ages by d/c. 



From Fig. [HI it clearly appears that all sources do not 
contribute the same way. We first list here only the most 
important ones: 



ay = 0.70 ± 0.10, and az = 0.81 ± 0.16. They also gave 
brightnesses at 1 GHz of each parts: Sx = 1, 160 ± 200 Jy, 
Sy = 440 ± 240 Jy, and Sz = 400 ± 260 Jy, which imply 
that SxYZ = 2, O OP ± 700 Jy c onsistent with the value of 
1,750 Jy given bv lGreenI (l2009l) . 



G203.0+12.0-3.3 Also called Monogem ring, this object is 
considered as a "probable remnant " by the Green catalog . 
Results based on observations from lPlucinskv et al.l (|1996D 
derived a distance of ~ 300 pc an d an age of 86 kyr. Parallax 
measurements by iBrisken et al.l (pOOSl for the associated 
pulsar gave a more precise distance of 288]'^27 P*^- Being 
25° in diameter in the sky, it is possible to infer neither a 
correct power spectrum nor a brightness. 

G263.9-3.3 This SNR is one of the brightest radio object 
in the sky. Its large apparent size was the reason it was 
long considered as three different objects of the Vela con- 
stellat ion, and it is often called Vela (XYZ). ICha et ahl 
(jl999t) estimated its distance to be 250 ± 30 pc from ab- 
sorption measurements. This is consistent with the paral- 
lax measurement of the associated pulsar (PSR) B0833-45 
from ICaraveo et al.l (|2001[) which gave 294l5p pe. Pulsar 
period derivative a nd nebula evolutio n implied the sam e 



age of - 11,200 yr (|Tavlor et al 1993i: iMiceh et al.ll2008[) - 
Radio observations reviewed bv lAlvarez et al. ( 2001[ ) gave 
the three spectral indices of the structure: ax = 0.39±0.03, 



Antlia Although not yet in the Green catalog, this remnant 
w as discovered by McCuUough et all ()2002[ ) and confirmed 
bv lShinn et~all (I2007D . It lies at coordinates 276.5+19.0 and 
at a distance of only 60-340 pc from the Sun (Sedov estima- 
tion). Other considerations (maximal size, column density 
and ^^Ti abundance) favors the lower bound of the dis- 
tance. Its age is probably more than 1 Myr. However its 
angular diameter is so large (24°) that a correct spectral 
index cannot be inferred. 

Hereafter we list other nearby SNRs (<2. kpc). Some 
authors did not correct their age estimations by the light 
travel length, sometimes leading to results outside our light- 
cone. Therefore, a supplemental uncertainty of Ad/c in the 
age should be added in some cases. 



G18. 9 5-1.1 Using HI absorption observations. iFurst et al.l 
(|1989[ ) estimated its distance to the Earth to be 2 kpc, 
but they could not exclude the possibility that it could be 
an extremely and unu sually bright o bject located 15 kpc 
away. The study from iHarrus et al.l (|2004[) estimated the 
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age of the object at 4,400 to 6,100 yr (without taking into 
account the distance uncertainty). However, this value is in- 
consistent with causality constraints, so wc added the value 
correspon ding to the d i stanc e 6.5 kyr. Using ROSAT a ob- 
servation, [^ersFeOD (|1997f ) found a spectral index of 0.28 
but the remnant is hosts a great deal of substructures that 
makes this value difficult to appreciate. They also infer a 
luminosity of order 40 Jy at 1 GHz. 

G65.3+5.7 Sometimes this remnant is also called 
G65.2-I-5.7. The radio power spectrum is 0.58 ± 0.07 
from 83 MHz up to 4.8 GHz (IXiao et al.ll2009D . and the 
estimated distance varies from 800 pc to 1 kpc. In earlier 
studies, in the conte xt of a pulsar survey in several SNRs 
(jGorham et al.lll996[) . among which G65.2-I-5.7, the flux is 
estimated to 52 Jy at 1 GHz, but this l ast quantity is rough 
estima ted. Optical observation from iMavromatakis et al.l 
(|2002f) leads to an age of - 26 kyr. 

G67. 5+1.2 Also called D A 495, this object is highly polar- 
ized. iKothes et al.l ()2006D suggested that point-like sources 
may alter the measurement of the power spectrum, leading 
in that case to a lower limit of 0.38 ± 0.08. The bright- 
ness at 1 GHz they obtained is 4.6 ± 0.2 Jy, which should 
be considered as an upper limit since the re mnant seems 
to sit on a diffuse emission plateau. However. iKothes et al.l 
(|2008f) proposed that it is not a shell-type SNR but a PWN. 
Indeed the spectrum has a power break at ~ 1.3 GHz and 
no shell. The power spectrum is 0.45 ±0.1 and the spec- 
tral break is of 0.42 ± 0.22. The corresponding brightness 
is of order ~ 5 Jy at 1 GHz. They also give a distance of 
1.0 ± 0.4 kpc and an age of 20,000 yr. However, if there is 
a pulsar in this remnant, it has not yet been detected. 



G69.0-h2.7 More often called CTB 80, this remnant is es- 
timated bv lCastelletti et al.l (l2003l) to be ^ 20,000 yr old 
Multi-frequency review from iGastelletti fc Dubne^ (|2005[ ) 
gave a spectral index of 0.36 ± 0.02 and a brightness of 
65 J y at 1 GHz. Howeve r, the Canadian Galactic plane 
survey (|Kothes et al.ll2006^ gave a spectrum of 0.2 ± 0.1 
and a brightness of 60 ± 10 Jy. Many authors cite 
iStrom fc StappeTi (l200l as a reference for a distance of 
^ 2 kpc based on H I absorption observations. The study by 
IZeiger et al.l (l2008f) seems to reveal an association with the 



pulsar B1951-I-32. However it might be much older (51 kyr) 
than the remnant. 



G74.0-8.5 Also known as Cygnus Loop this ~ 10,00 yr 
old remnant has been shown by iBlair et "aD (|2009[ ) to 
be farther than 576 ± 61 pc away. This is consistent 
with H ubble Space Telescope obse rvations fro m Blair et al] 
1200^, which found 540l8^° pc. ISun et all (I2006D found 
a spectral index of 0.40 ± 0.06 and an integrated flux of 
175 ± 30 Jy at 1 GHz. 



G78.2-f2.1 This remnant is also called 7 Cygni. The 
catalog based on the CGPS gives a brightness of 
27 5 lb 25 Jy at 1 GHz takin g into account the work 
of iLadouceur fc Pineau E (|2008f ). which yields a very pre- 



by iMavromatakid (|2003[ ) revealed an age of ^ 7,000 yr. 
However, it is not possible from HI observation to infer a 
distance, depending on the method they find either 1 or 
4 kpc with large error bars. The value of 1.5 kpc s eems 
to preferred by most authors (see iBvkov et al.l |2004 and 
references therein). 



G8 2.2-I-5.3 Also called W63, this remnant is estimated 
by IMavromatakis et al.l ()2004l ) to be at a distance be- 
tween 1.6 and 3.3 kpc and an age between 13.5 and 26.7 
kyr. However, it is not clear whether the Sedov analy- 
sis they h ave performed is lici t in th is case. Previous re- 
view from iRosado fc Gonzalei (|l98lh suggests a distance 
of 1.6 ± 0.3 kpc. The catalog based on the CGPS gives a 
brightness of 105 ± 10 Jy at 1 GHz and a spectral index of 
0.36 ± 0.08. 



G89.0+4.7 According to lBvun etld] ()2006[) . this remnant, 
also known as HB 21, is at a distance of 1. 7 ± 0.5 kpc from 
the Sun. The review bv lReich et all (|2003D suggests a spec- 
tral index of 0.41 ± 0.02 and a flux density of 228 ± 5 Jy 
at 1 GHz. However, the more recent CGPS gives an in- 
dex of 0.27 ± 0.07 and a br i ghtne ss of 200 ± 15 Jy. A 
study from iLazendic fc Sland ()2006() estimates the age to 
be 5,600 ± 280 yr but they used the for mer distance es- 
timation of 0.8 kpc by iHumphrevd ([lOT^) ; it is not clear 
how this age would change for a larger distance of 1.7 kpc. 



G93.7-0.2 Also referred to as CTB 104A or DA 551, 
this highly pol a rized and large object was studied 
by IKothes et al.l (|2006D . They measured a spectral in- 
dex of 0. 52 zfc 0.12 and a b rightness at 1 GHz of 
42 ± 7 Jy. lUvaniker et al.l (I2002D es t imate d its distance to 
be 1.5 ± 0.2 kpc but in lKothes et all ()2006l ). citing the same 
reference, the authors give 1.4 k pc. The only age estima - 
tion we could find is the one from iMantovani et al.l (jl982l ). 
which suggests 29,000 to 74,000 yr. However most of the 
other results of these observations have been corrected by 
later works and therefore it is not clear to us how trustwor- 
thy this age estimation is. 



G114 .3-I-0.3 Although some authors (see e.g. iKaolan et al.l 
[2006h stih use th e large value of 3.5 kpc, the study of a 
HI association by lYar-Uvaniker et al.l (|2004D gives 0.7 kpc. 
Same authors also esti mated an age o f 7,700 yr. Because of 
its proximity to Gas A, IKothes et al.l (|2006f ) were not able 
to measure the flux and the spectral index. However, they 
give, as an average of previous works, a spectral index of 
0.49 ± 0.25 and a flux of 6.4 ± 1.4 Jy at 1 GHz. 



else spectral index of 0.75 ± 0.03 correctly taking into 
account the thermal subtraction. A deep optical survey 



G116 . 5-1-1.1 While some authors (see e.g. iKaplan et al.l 
l2006t ) still use the lar ge value of 4 kp c, the study 
of a HI association by lYar-Uvaniker et al] (|20b4il gives 
1.6 kpc. Same authors also estimated its age to be be- 
tween 15,000 and 50,000 yr, they wrote that they tend 
to trust mo re the young e r age . We decided to assume 
20 ± 5 kyr. iKothes et al.l (|2006t ) measured a spectral in- 
dex of 0.16 ± 0.11 and a flux of 10.9 ± 1.2 Jy at 1 GHz. It 
is possible that the pulsar B2334-I-61 is associated with this 
remnant (or to the next one), although the age and distance 
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estim ations given by the ATNF catalog ([Manchester et aH 
1200511 are not in close agreement. 



Gl 16.9+0.2 Also called CTB 1, this remnant is consid- 
ered by some authors (see e.g. iKaplan et al.l |200"6[) to lie 
3.1 kpc away. However t he study of HI association by 
lYar-Uvaniker et al.l (|2004 ) gives 1.6 kpc. It seems that it 
is very close to G116. 5+1.1. The same authors also esti- 
mated its age to be between 15,00 and 50,000 yr favoring 
a younger age. lKothes et all (|2006( ) measured a spectral in- 
dex of 0.33 ± 0.13 and a flux of 7.9 ± 1.3 Jy at 1 GHz. 



Gl 19.5+10.2 Also known a s CTAl, th i s rem nant is the 
first object seen by Fermi (lAbdo et al.l l2008l ) . The most 
recent estimat i ons o f its age and distance are those of 
iPineault et al.l ()1993() . who estimated the distance to be 
1.4 ± 0.3 kpc and th e age to be between 5 ,000 and 15,000 
yr. The later study of lPineault et al.l (|l997t) revealed a spec- 
tral index of 0.57 ± 0.06 and a brightness of 40 ^ 45 Jy at 
1 GHz. 



age of the pulsar. I Reich et"al] ([2003) measured the bright- 
ness of the sourc e at 863 M Hz. Using the spectral index of 
0.30+0.15 fromH iao et al., (.2008i) , one finds a brightness 
of 74 ± 12 Jy at 1 GHz. lXiao et al.l (|2008f ) also showed that 
the spectrum is broken at ~ 1.5 GHz above which the index 
increases to 1.35 ± 0.20; this is probably due to the diffuse 
component. 



G184.6-5.8 One of the most famous, more known as Crab 
nebula or 3C144, this is the remn ant of the h is torica l 
Supernova SN1054. As explained in iKaplan et al.l 



even though this object is used to calibrate many instru- 
ments and hence very well studied, its exact distance can- 
not be measured precisely because its extreme brightness 
prevents from a parallax study. The nominal distance in 
the literature is 2.0 ± 0.5 kpc. The corresponding age is 
therefore 6,000 to 9,000 yr. According to the Green cata- 
log, its spectral ind ex is ~ 0.3 and i t brightness is 1,040 Jy 
at 1 GHz. However . I VinvaikinI ()2007l) showed that this value 
is decreasing with time. 



G 1 27. 1+0.5 Also known as R 5, this remnant was studied 
bv lLeahv fc TiarJ (|2006[ ). who used radio observations from 
the CGPS. They determined an age of 2-3 10^ yr, a spectral 
index of 0.43 ± 0.10, and a brightness of 12 ± 1 Jy at 
1 GHz. For the distance, thanks to HI line survey, they 
found 1.15 kpc which is consistent with an association with 
NGC 559, located 0.9-1.13 kpc away from the Sun. 



G156.2+5.7 Disco v ered during ROSAT survey by 
iPfeffermann et aTl (|1991l) . this remnant is one of the 
faintest in the G alaxy with a brigh tness of 4.2 ± 0.1 Jy at 



Reich et al. 19921 ). In t heir study of the 
Pannuti k AllenI (|2004[) found that the 



1 GHz (see e.g. 
X-ray emission, 
photon emission is in agreement with an electron density 
with a power law of 2.0 j^J'y and a cut-off at about 10 
TeV. Its distance from the Sun was first estimated to 
be a few pc, although it seems from optical study by 
iGerardv fc FesenI ()2007() that the remnant may be as close 
as 300 pc and younger th an 10,000 yr. In their latest 
paper. iKatsuda et al.l (|2009[ ) use 1 kpc. 



G160.9+2.6 Also called HB9, iLeahv fc TianI (|2007t ) using 
the CGPS and HI observations estimated a spectral index 
of 0.48 ± 0.03, and a distance of 0.8 ± 0.4 kpc. Concerning 
the age, they conclude that the Sedov age for HB9 is 6,600 
yr and the evaporative cloud model yields ages of 4,000 
- 7,000 yr. However, the age of the corresponding pulsar 
B0458+46 seems to be 7 kyr. The brightness is of the order 
of 75 Jy at 1 GHz. 



G180.0-1.7 This bright radio remna nt is also c alled S147. 
The 143 ms pulsar PSRJ0538+2817 (|Andcrson e t al. 1991) 
is located within S147 and believed to be associated as the 
distance estimates for both the remnant (0.8-1.6 kpc) and 
the pulsar (1.2 kpc) agree. The characteristic age of the 
pulsar of 600 kyr is much o lder than th e estimated 100 
kyr based on radio data from ISofue et al.l (0,980.1 However 
the complexity of the environment of this remnant makes 
any age analysis extremely difficult. We therefore adopt the 



G189. 1+3.0 Also named IC 443, this 30,000 yr old 
remnant i s well-known for it s rich chemical composi- 
tion (see Neufeld et al.l I2007L and references therein). 
IWelsh fc Saih^ (120031) have calculated a distance of ^ 
1.5 kpc. The review bv lMufson et al.l (I1986D gives a power 
spectrum of 0.36 ± 0.04 and iReich et al.l (|2003D provide 
an intensity of 16 ± 5 Jy. An XMM-Newton observation 
(jTroia et al.ll2008[) implies a much younger age of 4,000 yr. 



G205.5+0.5 Also known as Monoccros Nebula or the 
Monoccros Loop, this S NR is believed t o be associated with 
the Rosette Nebula (see lOdegardlll986D which lies at a dis - 
tance of 1.6 ± 0.2 kpc from the Sun (|Bonatto fc Bical[2009l) . 
This special position makes it a very interesting probe 



for cosmic ray acce leration study (see iFiasso n et al. 20081). 
iLeahv et al.l (I1986D have estimated the age of the object to 
~ 29,000 yr. However, this value is ver y mo del- dependent , 
and ultraviolet spectroscopy analysis bv lWelsh et al.l ()2001|) 
suggests t hat it is even older (up to 150,000 yr). The 
review by iGraham et al.l (|l982l) suggests a spectral in- 
dex of 0.47 ± 0.06 for the radio emission and a flux 
density of 156.1 ± 19.9 Jy at 1 GHz. The results of 
iBorka Jovanovic fc UrosevicI (|2009( ) which have higher sta- 
tistical significance prefer an index of 0.66 ± 0.20 with a 
brightness at 1 GHz in agreement with previous ones and 
a distance of 1,630 ± 250 pc. 



G266.2-1.2 Often referred to as RX J0852.0-4622 or Vela 
Junior, the detection of the radioactive decay line of ''^Ti 
seems to prove that it is a very young object. However its 
proximity to Vela makes any estimation of it brightness and 
spectral index quite difficult. Based on an evolution study, 
iKatsuda et al.l (|2008l) estimated its age to be between 2,700 
and 4,300 yr and its distance to be ~ 750 pc. It is a little 
puzzling that no historical record of this supernova explo- 
sion, that happ ened around year 1, 300, has been found. It 
was proposed in lRedman fc Meaburn (2005) an association 
with pulsar PSR J0855-4644, based on a reestimatc of the 
pulsar distance. 
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G315. 1+2.7 The discovery by IStuEareFaD (|2007D of this 
extremely large remnant revealed a distanc e of ^1.7 kpc. 
Reana lyzing the former radio observation by iDuncan et alj 
(jlQOTT) they found a spectral index of 0.7. However, it seems 
that some ambiguity remains concerning the brightness. 
The age is not given precisely but it is considered as old. 
Without anymore precision we have taken 50 ± 10 kyr. 



G330. 0+15.0 Also call ed Lupus Loop, this remnant is 800 
pc away from the Sun (iKaplan et al.ll2004|) . However, a re- 
view bvlShinn et all (l2006l^ gives 1.2 ± 0.4 kpc. According 
to iKaplan et al.l ( 20061 ) (and references therein) it is 50 
kyr old. Very few data are available for Lupus Loop in 
radio wavelength. To our knowledge, the latest arc from 
iLeahv et al.l (|l99lh , who concluded that its spectrum could 
not be reproduced by a single power law. 



G347.3-0.5 Associated with the gamma source RX J1713.7- 
3946, this SNR lies about 1 kpc away. The age of the 
first light on Earth is estimated to be ~ 1,600 yr, con- 
sistent with the historical Chinese record of an SN ex- 
ploded in AD 393 (see e.g. iMorlino et alll2009t) . For a dis- 
tance of ^ Ikpc away, this me ans an age of ~ 4.9 kyr. 
in ison fc Cassam-Chenail (j2005f ) explained that because its 
environment is extremely complex, it is very difficult to sub- 
tract the background and to fit both a power spectrum and 
the brightness. 



